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ABSTRACT 


= A-displacement solution for the calculation of structural influence 
coefficients (SIC's) is presented, The formulation utilizes the lumped 
parameter approach that is consistent with collocation flutter solutions. 
The structure is synthesized as concentrated mass elements connected by 
massless elastic plates and/or beams. There are two methods of generating 
the mass matrix; they are: 1) lumped concentrated mass points, 2) consistent 
mass matrices, Along with the calculation of the SIC's, the natural 
vibration modes and frequencies are calculated, There are two options 
for punching out the flexibility matrix for use in subsequent COFA computer 
programs. Option 1, punches out the full flexibility matrix; Option 2, 
punches out the reduced flexibility matrix eliminating the rows and columns 


pertaining to structural attach points, | 
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1.0 INTRODUCTION 


In order to determine the aercelastic behavior of a wing or control surface, 


it is necessary to know the aerodynamics, elastic properties and mass distribu- : 
tions of the structure. The overall aeroelastic analysis is usually divided: into four: 
separate parta as shown in Figure 1. 4 


owe be 


Mass q Stiffness 
Distribution Distribution 


Aerodynamics 


Flutter 
Speéds 


Figure 1. Analysis Procedure 


This portion of the report describes the computation of the mass and 
stiffness distribution. The geometry of a wing or tail surface is too complex 
for the successful use of closed form analytical techniques. Therefore, a 
numerical type of analysis must be used. The end product of this analysis is 
the generation of overall influence coefficient and mass matrices referred to 
a set of node points arbitrarily picked on the surface of the structure. The 
finite element method (see Refs. 2 and 3) was used to form the required matrices 
for a planar structure. This technique is especially suited to solve complex 
structures and as used in the analysis ie general enough to handle the following: 

1. Combinations of beam and plate elements 

2. Arbitrary boundary conditions 

3. Lumped or distributed stiffnesses and masses 
A discussion of the theory and computer proyram which calculates the influence 
coefficient and the mass matrix as well as the structural modes and frequencies 


is given in the following sections. 


2.0 NOMENCLATURE 


C = Unknown Boundary Constants 

D = Plate Ridigity Constant 

E = . Modulus of Elasticity t 
F = Force ) 


K = Stiffness Coefficients 


M = Bending/Torsional Moment 
p = Pressure 
T = Coordinate Transformation 
t = Thickness 
w = Linear Displacement in x direction 
X, Vy 2 = Coordinate Axes 
s = Linear Displacement 
pa = Curvature 
e = Density 
T = Stress 
Y= Poisson's Ration 
& = Partial Derivative - 
[ ] = Square Matrix 
{ t = Column Matrix 
iL J = Row Matrix 


j 
f 3.0 TECHNICAL DISCUSSION 
3.1 Influence Coefficients 
| The stiffness method approach is firse used to obtain an overall stiffness 
matrix of the structure. This matrix is reduced by partitioning and then in- 
verted to obtain the influence coefficients at any. desired set of control points. 
: The number of control points are denoted by N. At each node, three degrees of 
freedom are specified: two rotations and the normal displacement. Therefore, 
a stiffness matrix of approximately 3N degrees of freedom is first formed by 
superimposing individual plate and plane grid beam element global coordinate 
matrices. The matrix will be somewhat smaller than 3N degrees of freedom since 
boundary restraint conditions will reduce the size of the matrix. To illustrate 
the matrix condensation method used in the computer program, we will assume that 
we have N control point normal displacements and M displacements which must be 


eliminated. The overall stiffness matrix is given as 
iK] = (1) 


The structural equilibrium matrix equstion can be written as 


[ K K § F 
NN NM N : N (2) 
Pai Kis bu 4 
3 


We now assume that forces at the points to be eliminated are small and can 


be neglected. Therefore, 


K WN K NM by a 
Kuan Kaw | 1 on : 0 
+ Hull = 
lel + P83 = £9 
{54} = ~[Kvna] [Kn] {0} 


(ul ~ [Xml al” [uel) {6 


and since - 
Mun] = [vn 


we have ' 


‘ £5} = (and = [Kun [Anne] [Aun]) {Fa 


If we now let T oe y 


[Fond = ( [Kal ~ [Kon 


then Equation (4) can be written as 


Varad 
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The matrix fig is called the structural influence coefficient matrix, 
The application of loads at the control points yield displacements at the 
control points by carrying out the matrix multiplication indicated in 
Equation (5). ' 
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Figure 2. Plane Grid Beam Local and 
Global Coordinate System 
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The computer program FLUENC carries out the required operations to obtain 
the influence coefficient matrix ES . A detailed description of the program 
can be found in Section 4.0. The program is written to form a 50 x 50 influence 
coefficient matrix. The influence coefficient matrix is punched out on cards in 
a format compatible with the Collocation Flutter Program. 

The plane grid beam global coordinate stiffness matrix used in the program 
was obtained from Reference 1 and is given in Table 1. ‘The local and global 
coordinate systems are shown in Figure 2. The figure also contains the sign 
convention for the six dezrees of freedom for each element. 

The triangular plate stiffness matrix given in Reference 2 was used in the 
computer program. The plate element can be materially or geometrically ortho- 
tropic as treated in Reference 3. Stiffened' plates can be considered to be 
geometrically orthotropic. ‘The sign convention and nodal degrees of freedom are 


shown in Figure 3. 


z' parallel to z 


x', y' - follows principal 


axes of material 


- angle between x' and 
x axes 


Figure 3. Orthotropic Triangular Element 
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Following the analysis given in Reference 2, the deflection shape of the 
plate element is assumed to be of the form 
3 t 
we C4 C,H + C39 4+ CE Ks Coxe er Coy 
2 % x 3 
+ C,x?+ Co (ay + xy) + Cyy 


or (6) 


ur = [N] | Cc i 
The unknown constants Cc» Ci, °°"; Cy can be written in terms of the nodal 
cisplacements by 5, ---, b, by using the boundary conditions 
ur = 8, 
at X20, y= C Gur / Og = $. 
der | In = - 3; 
UP = Da 
at KO, oe ) Ie] 24 = Se (7) 
| Qeur/ Ox = 5, 
ur = by 
at HRs, 9° 9s der / dy + bs 
de [ox ze 53 


Using Equation (6) in conjunction with the boundary conditions given by 
Equation (7) yields 
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The constant vector {c } can be obtained in terms of the nodal displace- 
ments by inverting the matrix [c] . Therefore, 


fe}= [elf » 


The curvatures for a Flat plate element are given by 


eS 


Ex ae | de® 
ie} = 6, aie Dur / Ay - 
exy 2 Our / dx Oy 


Substituting Equation (6) into Equation (10) yields 
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where i 
0 C 0 -2.0 0 -te  -2T ) 
CO ee Oo OO. Qn 70o 82 SO aoe ocho 
0 O 0 C -2 0 oO (4X44 O | 
Substituting Equation (9) into Equation (11) yields i 
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1f initial strains are neglected: then the moisent-curvature relation- 
ships can be written in the fora 


fo}= (mM, }= [oife} (3) 


where 


ID] =| Lb DB 0 (14) 
7) 0 Dy 


for a materially or geometrically orthotropic place. For an iaccropic plate 
Equation (14) reduces to 


? / v O - 
Ss bee he 0 (15) 
te(1- v*) O l-2 
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The [p | matrix must undergo a transformation if the principal axes of 
the material do not coincide with the local coordinate axes. The components 
of strain in one coordinate axes system are related to the components of strain 


in anether coordinate axes system by the matrix equation 


fe} = [7 Tf e] a6 


(The prime refers to the components of strain referred 
to the x'-y' axes in Figure 3) 


where 
r Cos *B omy'7? . 2 tunyse “of3 
i ] 2 ye te Zesryptcosfd 
LT = Sin COS 2 2 / /- (17) 
Py , 4 4g 
sin f3 cas 8 - genss cor f3 i) a Ed 


Since the internal work is constant no matter which coordinate system 


fe} fe} - fo} fe} 


is used 


or by Equation (13) 
fe") Lo’ fe] 


and by using Equation (16) 


feVEr][o'1[ 7] fe} = fe} [2] {6} 


Therefore 
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The stiffness matrix for a typical element (1) GB) is given by 
T 
[K] = f[lallollal wy 
A 


or by Equations (12) and (19) 


[k]- (eT (floF U7 foUtr} Te] arey)bel 


Now let 


(5) = fffol iro llrl Te] wy 


and carrying out the indicated matrix multiplications ylelds 


[D | Jf (se 7able é) dydy 
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(21) 


(22) 


In order to simplify the integration required for evaluating the matrix 


in Equation (22), it is suggested in Reference 2 that the independent variables 
be changed as shown in Figure 4. . 
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Figure 4 
Coordinate Transformation 


The relationships 


E(l-a)4, 
El(1-7)9 ny, | (23) 
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are used for the change of variables. The terms in Equation (22) can’ now 


be evaluated by using the relationship 


I(x", 9’) = ff xy" dye dy 


Tia, 9") « ff wy" |Tlx,y)| dg fe (24) 
where 


J (x, 4) = 


(25) 
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Substituting Equation (23) into Equation (25) yields 
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T(x,4) = § As Ye 


Substituting Equations (23) ane (26) into Equation (24) yields 


T(x", y") = j f pena (a) 4 
ty) as" 4, ds dy 


which can easily be evaluated for any m and n. 
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3.2 Mass Matrix 


D'Alenbert's principle can be used for the formulation of the mass 


matrix. If masses are attached to the nodes of the structure, then the 
nodal dynamic forces are 


us) 


[P}-- 
& 
dé 
: where 
M, 0 
M, 
[mM] = rs 
0 M, 
is a diagonal matrix. The mass of beam and plate elements are usually dis- | : 
tributed throughout the structure. Therefore, the distributed pressure bas | 
loading can be written in the form ] 
t y 
p--p Ze oof 
d¢ 


Substituting Equations (6) and (9) into Equation (30) ylelds 
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Since the equivalent element nodal forces can be computed from 


the equation 
é g 
{P}" = -Sfel'p dV oF 
thea 
[Py =| [lal'Tel ear] { 5} (33) 
Therefore the elemental consistent mass matrix is given by 


[m] - {[R] La] eav a 


The consistent mass matrices given in Reference 2 (see Tables 3 and 4) 
are used in the computer program. 

Once the elemental consistent mass and/or lumped mass matrices are 
computed, then the overall matrix is obtained by foilowing the same technique 
as used in assembling the overali stiffness matrix. 

The overall mass matrix is reduced by using Equation (3a). We again 
assume that we have N control point normal displacements and M displacements 


which must be eliminated. The overall mass matrix can be written in the form 


M M (35) 


and the displacements 


{5 | ae (36) 
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From Equation (3a) we have j 
{ 544 . Kane | Ky | 15y{ SI 


Since the virtual work of the reduced mass system must equal the virtual 
work of the true mass system 


- {45,97 [M,] {Sy} = - {45} Tm] {5} : 
where . 
144, = virtual displacements of control points 


{as { = virtual displacements of complete system 


[M, ] = overall reduced mass matrix 


Equation (37) can be rewritten in the form 


Se eT ee J 
(as, OMI {ayd = [aby as, fm] 4 ' 


Substituting Equation (3a) into Equation (38) yields ms 
(45,7 TMH St f45T [2 [Kelle EMI, oe Hal 
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which yields the result 
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The reduced mass matrix given by Equation (39) is calculated in the | 


computer program. 
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3.3 Modes and Frequencies 
Since the design engineer may find it useful to know the mode shapes 


and natural frequencies of the structure, this information can be obtained 


by using the NMODE option in the computer program. If no external forces 


toed 


. 
i 


are present then the reduced mass and influence coefficient matrices are 
related to one another by the relationship " 
1 

= t oe d ; 

Es 1S, - -[m]f 5) (40) | 

! 

For determining natural frequencies, the deflections 15,5 can be : 
written es 5 
ty 


[by } = 5} sin tot (41) : 


Substituting Equation (41) into Equation (40) yields 
=e 2 4 
el { &.} “ w [MJ {5} (42): 


The solution of Equation (42) yields the natural frequencies, ‘¢ , and 1 


the mode shapes {5.5 . Since fel and [u, | are both symmetrical 
matrices, the mass matrix M.| can be triangularized 


[MJ] = [2 II be (43) | 


where 


a ee ee ee ee 


Substituting Equation (43) into Equation (42) yields 
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Equation (44) may be written 


PTET okTsy on 


or 


[A]{B} = ot {8] as 


where [A ] e L ie [ fon) | uy 
(S.J -[e] [4] 


An eigenvalue subroutine using the Givens method was used in the 
computer program package to solve Equation (45). The Givens method is 
fully described in Reference 4. 


Note that the dynamical matrix [A} in the form described above is real 
and symmetric which is required by the Givens method. Conventently, 
fu) and (."] are in triangular form which is used in the computer program 


package to save core storage space. 
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4.0 PROGRAM DESCRIPTION 
H 

Computer program FLUENC written in FORTRAN IV carries out the operations 
set forth in Section 3.0 for generating the structural influence coefficients 
and mass matrices required by the Collocation Flutter Program. Briefly, the 
structure is assumed to be representable by a planar network of beams and. 
triangular plate elements connected at discrete joints. At each joint, if 
there are no restraints, the program assumes three degrees of freedom; that. 
is, one displacement normal to the plane of the structure and two rotations. 
The program first synthesizes the stiffness and mass matricee for the entire 
structure including all degrees of freedom from the data input for the beam and 
triangular plate elements and from the restraint information input for the 
joints, It then reduces the stiffness and mass matrices by eliminating all 
the rotational degrees of freedom and leaving only the normal dieplacements. 
As a final step, the program inverts the reduced stiffness matrix to obtain 
the influence coefficients. 

Other features of the program include the option to compute lumped masses 
or to compute the consistent mass matrices for the beam and triangular plate 
elements or both, Also, the triangular plate elements may have either 
isotropic or orthotropic properties, There is an additional option to expand 
the reduced frequency matrix to include the degrees of freedom representing 
the restraint joint (one joint on a movable surface; two joints on a fixed 
component). This is accomplished by adding one or two zero rows and columns 
to the reduced flexibility matrix corresponding to the mass numbers of the 
attach points involved. 

In the sections that follow detailed instructions are given for the 
preparation of input data and a description is given of the output illustrated 
with several sample problemeg, Aliso included are listings and flow charts of 
the program and a discussion of the processing requirements. 

4.1 Description of Program Input 

The following instructions describe the input data, their physical units, 
and the FORTRAN format they must be punched with. The input quantities! names, 
all in capitals, are their FORTRAN names and, for reference, their equivalent 
names in Section 3,0 are listed in Appendix D. 

4.1.1 Title Card, format (12A6) 


Two cards; any alphanumeric statement {n columns 1 to 72, 
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4.1.2 Problem Size and Control Information, format (715) 


colum | 1-5] 6-10] 11-15] 16-20 | 21 - 25 
NJTS = number of joints in structure (50 maximum) 
NR = number of joints with one or more restraints 
NBE = number of beam elements in structure i 
NPE = number of triangular plate elements in structure 
NM@DE = number of eigenvalues and eigenvectors desired (9 maximum) i 
MKEY = 1. do not compute consistent mass terms for beam and/or 
triangular plate elements 
= 2. compute consistent mass terms for beam and/or triangular i 
plate elements 
NLUMP = number of lumped masses input. Only lumped masses corres- / 
ponding to the normal displacement at each joint may be 
input. 
4.1.3 Material Properties | 


(a) Number of Materials, format (15) 


NMAT = number of materials for which properties are input (10 max.) 


(b) Properties, format (4E10.3) 


Input NMAT number of cards, one for each material. 


eye ame eanremen 


Column 11 - 20 
Name YM (1) DENS (1) 


erat os oon ce-enwe memoen! 


Young's modulus of elasticity divided by 109; pai 


YM(i) = 
PR(i) = Poisson's ratio 
GE(i) = modulus of rigidity; psi. If input as 0, it will be computed 
from the following formula: 
YM(i) 
c(i) = -—_—_—_—_—_— 
2 [1 + PR(1)] 

DENS(i) = material density; 1b/in?, Not requlred 1f MKEY = | 
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4.1.4 Joint Coordinate Cards, format (LOX, 2E10.3) 


Input NJTS number of cards, one for each joint. Also, the structure 


is assumed to lte in the x-y plane. 


i 
l 


m = joint number (must be input consecutively starting with 1). 
May be placed anywhere between columns 1 and 10 


X(m) = x coordinate of joint m; inches 


Y(m) = y coordinate of joint m; inches 


4.1.5 Joint Restraint Information, format (415) 


Moths sth6to se Sei ae 


Input NR number of cards, one for each joint with one or more restraints. 


Column 

Name 

JT = number of joint having one or more restraints 
Ml = O free in the z direction 


= 1 fixed in the z direction 

M2 = O free to rotate about the x axis 
= 1 fixed about the x axis 

M3 = O free to rotate about the y axis 
= 1 fixed about the y axis 


i Ci oacti aah ips ale eich ta ANGEL DS CARES Piya BEAT Yee 18 Stet Sh 


4.1.6 Lumped Masses, format (15, 5X, E10.3) 


? 
Input NLUMP number of cards, one for each lumped mass. 


E 
E 
i 
i 
| 
[ 
| 
| 
| 


Column 

Name 

JMASS = number of joint for which lumped mass La input 
= lumped mass, Ib. 


| RSMASS 


If more than one lumped mass is input for a particular joint, the 


| program will sum the masses. 
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4.1.7 Beam Element Properties, format (3E10.3, 315) 


Input. NBE number of eoree one for each beam element. 


~ Pepe pe] = [me] | 


AR * area of beam cross section, in? 


XI ™ moment of inertia of area, in’ 


YJ = effective torsional moment of inertia, in’ 


MAT 


material code corresponding to one of the materials input 
under paragraph 4.1.3. 


JTNR, JIFR = joint numbers at the ends of the beam element 


4.1.8 Zrigngular Plate Element Properties, format (£10.3, 515) 
Input NPE number of cards, one for each triangular plate element. 


PTH = plate thickness, in. 


MAT = material code corresponding to one of the materials input 
under paragraph 4.1.3 


JT1, JT2, JT3 = joint numbers at the three corners of the triangular 


plate 
Restrictions: 
a) The order of the joint numbers must be given in a clockwise 
manner as follows: j7T3 
J yi, 
TT 
a 
— eZ JT2 


b) The angle fermed by the edges of the triangular plate at 
JT1 must not be 90". 


NDX = © the plate has isotropic propertles and the flexural rigidity 
terms are computed from 


YM(MAT) x PTH? 


pero 12 fi - [prowr) ] ?} 


Di = | PRAT) J x DX 
20 


——- myer (CC (tC 


os 


a 


y = 1 the plate has orthocrepic properties and the flexural rigidity 
terms are input by the next card [ format (4E10. 3)} 


l1- 201 21-30 | 31 - 40 | 41 - 50 | 
pe pape] 


DX, DY, Dl, DXY = flexural rigidity terms, in.1Lb. 


BETA = angle between material principal axes and the triangular plate 
locel coordinates as shown below 


YY 


Material Principal 
Axes 


Le X 


Triangular Plate Local 
\ iN Coordinates 


4.1.9 Option to Expand Reduced Flexibility Matrix 


Note: The following card (NCOD) is always required at the end of all input 


data for any one particular case, whether or NOT the option is to be executed 
FORMAT (116) 


Column 1-6 
Name NCOD 
item (1) 


NCOD = 0 Option not executed 


=z 1 Option executed 


IE NCOD = 1, the following card is required | 


FORMAT (318) 


Column \ i 
Name 5 
Item | 
NR = Number of boundary points used (1 or 2) 4 
NNE = Mass number of first attach point 
| NWO = Mass number of second attach point, if NR = 2 
NWO = O or left blank if NR = 1 


. We 
Perens 


i 2 


A Pe eaae See-t aaa tl 


To input more than one problem, the user need enly repeat the cards in 
paragraph 4.1.1 through 4.1.8 for each additional problem. 


4.2 Description of Program Output 

The program prints out all the input data for every problem followed by 

the solution consisting of the reduced upper right triangular stiffness 
{lb/in), flexibility (in/1b) and weight (1b) matrices as well as the modes and 
frequencies when these are requested on the card in paragraph 4.1.2. ‘The 
stiffness, flexibility, and mass matrices that are printed/punched out only 
contain terms that are associated with the normal displacement "z". This 

is done so that when the flexibility matrix is used in subsequent collocation 
flutter analyses only the easential degrees of freedom are included in the 
flutter analyses, Also, the matrices are reduced to eliminate control points 
associated with fixed points (boundaries). If it is desirable to include the 
boundary points, it is only necessary to intersperse rows and columns of zero's 
at the proper place in the matrices... Immediately following the joint restraint 
information in the output, the program prints out the coordinate numbers assigned 
by the program to the normal displacements at each unrestrained joint. The 
elements in all the reduced output matrices are ordered according to these 
coordinate numbers, 

In addition, the program punches out the entire flexibility and weight . 
matrices row by row with the format (1P6E12.5) which is compatible with the input 
requirements of the Collocation Flutter Program. Each punched matrix is identified 
by a little card as the first card. 


4.3 Sample Problems 


To illustrate the use of program FLUENCE, three sample problems are in- 
cluded in Appendix A. Each sample problem starts with a problem statement and is 
followed by a listing of the input data and the output of the program, The Firat 
sample problem is a simply supported uniform beam composed of five beam segmentn, 
The second is a uniform cantilever plate divided into 72 triangular plate 
elements, and the third is a lumped mass and beam network simulating a misulle 


control surface, 
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4.4 Processing Requirements 
Program FLUENCE has been run on the GE-635 computer and it required about 


31,000 cells of core storage. It is expected that the program storage require- 

ment will be about the same on other digital computers. In addition to using 

the input and output files, 05 and 06, which are standard for the GE-635 computer, 

the program requires six other peripheral files, five of which are designated in 

the program by the numeric codes 07, 08, 19, 10 and 11, and the sixth is the card 

punch file. ) , 
There is no general formula for determining the run time required for a problem, 

but if a GE-635 computer is used, an estimate may be made from the times required 

for the three sample problems in Appendix A, which are as follows: 


“No. No. of No. of Consistent | Lumped | No. of 
of Beam Plate Masses Masses jMcdes & Freqs 
Joints | Elements | Elements | Computed | Input} — Computed 


4.5 Program Listing and Flow Chart 
In the event future changes are needed in the program, a listing of the program 


is included in Appendix B. The program consists of a MAIN deck, 24 subroutines and 
one function subprogram. MAIN has the function of reading in data, numbering the 
coordinates (subroutine C@@RDN), generating the codes for assembling the stiffness 
and weight matrices and calling the subroutines which develop the stiffness and mass 
terms for the beam and triangular plate elements. When the entire stiffness and 
weight matrices have been established for the whole structure, the MAIN program calls 
a subroutine which reduces these matrices as discussed before and determines the 
modes and frequencies as well. 

The 24 subroutines and one function subprogram can be divided conventently Into 
five groups according to their function. The first group consists of those routines 
that develop the beam stiffness terms; these are TRANS and BEAMK. ‘The second group 
consists of the routines which determine the beam mass terms; these are ‘IRANS and 
BEAMM. The third group develops the triangular plate stiffness terms and these are 
PLATEK, CMAT, MINV, DINMAT, MATMPY, DMAT, DBLINT and PLYMP. ‘The fourth group deter- 
mines the triangular plate mass terms and these consist of PLATEM, CMAT, MINV, DINMIIM, 
MATMPY, DBLINT and PLYMP. The fifth group of subroutines reduces the stiffness and 


23 


+? 
1 


SoS ASRS AIR It a OLEATE AMEN 


i 
Fi 
; 
: 
; 
i 


and mass matrices, finds the eigenvalues and eigenvectors and outputs the solution. 


This group is comprised of EIGEN, VIVID, ZR@MAK, ZR@MAM, SYMINV, EIGMAT, BIGMAT, 
LOOP1, LODP2, LOPP3 and LAPS. 

Since the program Listing is annotated extensively with comment statements, 
no further explanatory remarks are given here for the program. However, to faci- 
litate the understanding of the interrelationships among the many subroutines, a 
flow chart of the entire FLUENC program is included in Appendix C. 
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Identical to Sample Problem 1 with the addition of lumped mass 
input at joint 3 and 4, 
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Sagple Problem No. 3 


Missile Control Surface Model 
(Modeled with beam elements and lumped weights) 


Find first five natural. modes ani frequencies. 


Note: Joint 11 is restrained from rotating about y 


Lumped Masses 


Joint No. 


i 
2 
3 
4 
5 
6 
7 
8 
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- Sagan ae ES OL TL TEL ETE TOTO LE TERE EASTON A eo aera ey 


I og t : 
: | ies, : 
Beam Element Properties. : 
mi , | Moment-of-Inertia ~ Torsional : 
Member 7 Area = Constant — : 
i i- j of : | 
4 

0.0030 

0. 6070 

0. 0070 

0.9029 

E = 3x10’ psi 
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MESSILE CONTROL SURFACE MODEL WITH 79 JOINTS 
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& FORTRAN DECK 
CHAIN PROGRAM FLUENC-FOR GENERATING STIFFNESS,FLEXIBILITY ant MASS 
tt MATRICES FROM PLANE GRID BEAM AND TRIANG. PLATE ELEMENTS 
NIMENSTON TITLE (293, 9M(18)- PRCTOD GECIN) DENS(T OUD XCSN) Yeu), 
TNAT C5) NRO CSD NRI6 5) C5) HF C50) (00). NOSE(9),0CSiz), 
LSTHOA, GB) SMMC SAD PLIK (9,9) o PLTIM( 9.9) »SSTEC11325),SM(31525}, 
LRSMASSE1509, 4091575), VALUC9), TEMR(SU) BCI 9K) CCE), DUMSCINY). 
GFCUSA, TY), 191m4(59), JBASS(59) 
INT+GER OUT 
EQUIVALENCE CSSTF C1 ),SM(4),405))9,0STH(1-.1)-SMM(1,!D-PLTKE1,1). 
LPM ¢1,.1)) 
Nat FARMAT(12A6) 
10s) FORMAT (1615) 
JO? -FORMAT(AFI0.4) 
INH FERMATCTOX.PFIN, 8) 
1924 FORMATOCSEIR. 3.495) 
INHE FORMATIE NN. S.995) 
NBA FORMATOCES, 5X. 690. 9) 
KANN FOPMAT(IHE , EYAL SIX IAAA) 
SANT FORMATC///6HNITS =94.9X%,6H NR =95,9X26H NRE =1559X+6H NPE =19.5X, 
T7HNMONE =13,5X.5HMKEY =J$3-5X./HNLUMP ef.) : 
nn? FNRPRATC///75HM A TER EAL PROPER TIF S shenteaenneeas 
Ler er eeereeacenveeree/ /XHNN. YNUNGeS MODULUS POISSON RATIO 
1 MORULUS OF RIRENETY DENSTTY,19C/S F200 K FD se OKs THD VUX DEL 46%, 
L4X%.F12.5)) 


HN’ FARMAT(Os//349H,) O FE KT CA ARH I NA T F S/35HJGINE NG. x 
1 ENORDG. yY roorRo.) 

Sant CHRMATCESL7Y.F INL re OXF INL) 

SUS FORMATOS//6-HA 0 FT NT RES YTRATING COD E- eeeuteeeene 


Leeacsseeeasves/H7THIOINY WO. Z DISPLACEMENT ROTATION ABOUT X 
{ ROTATION aARnQUT Y) 
Sau? FORMAT( 15,114,319, 928) 
SHR FORMATC///75H8 F A M FLEMEANT PROPERT IE S tan 
Lesarsapaesseoncaseseceees /FFHELEMFNT NO. A ! 
1 d AAT JOINT 1 JOINT P) : 
HHND FARMATCT6,BX.F9O.4 AX FSG LIX FOCAL AX, TP OAK 17, 9X. 17) 
BALAN FAPMAT(///17°HT RET AN BULL AR PLATE FLEFMEN T 
1 PRO PER FSF LE § Feat wma eOS ae eePeP SEE RETHEHHRTREE HER EEHEEH EHH ED 
Teese /t°PHELEMENT NO, T MAT JQUNT 1 JOINT 2? JUOENT 
1% DX DY v1 NXy BETA) 
5M11 FORMATCTO,8 ©FR4 TKO EAs GK EPOX LP oRND IP AX EDL Se NK, F110 AX, 
LFI1.68, TK, F1 1. 5X-FHA.7) 
5020 FARKAT(///6°HC NURI ENATE MUMBERS FOR FACH 7 DISPLACEMENT AT EACH UN 
IRESTRAINED oOENT//5HJIOINT NO. COORD. NO.) 
§N21 FORMAT(C15,914) 
AIF FORMATO///PrHL OM POE OD WF TRH T §.2SHJOINT NA. WEIG 
VHT) 
SNe FORMAT(I1%,6 -F10.4) 
Wvsu 
C RISE ASSEGNMENTS 
TNs 
CALL FLGEOF CIN. TV) 
NUT=6 


ei 82 


nn eR a ye ae nee mere cme ee eK Mia am ce apne no 


iver 224, fee: amare. ne ote uke a - - + 
Se tee 
a 


MNISCs7 
NNYSC=A 
THTSC=9 
JUTSE=10 ? 
KNTSC=1 4, , 
c REGIN INPUT OF DATA 
tee READCEN,TG°U) CTITLECT), 121,74) 
TECHVNE.U) CALL EXET 
REWIN) MDISC 
REWIND NDISU 
REWIND IDISC 
REWIND JBISC 
REWIND KDISv 
WRETECOUT,-naa) (TETLECE), $29,524) ; 
READCIN, 10%) ) NJITS. NR,NBE NPE, NMOUE> MKEY, NLUMP 
C NJTS=NO. OF JOINTS, NR=NO. OF JOINTS WITH RESTRAINTS 
c NHF=NG. OF KBFAM ELEMENTS, NPE=NO. OF TREANGULAR PLATE ELEMFNIS 
C: NMOLE=NO. OF EIGENVALUES AND EIGENVECTORS NESIRED 
c MKEY = 1 DO NOT COMPUTE ELEMENTAL CONSISTENT MASS TERMS 
c MKFY = ¢ COMPUTE ELEMENTAL CONSISTENT MASS TERMS. 
c NLUMP = NO. OF LUMPED MASSES INPUT 
WRITECOUT, 9007) NJITS, WRANBE.NPE.NMODE,MKEY, NLUMP 
f INPUT 4ATERIAL PROPERTIES 
READCIN,10eL) NMAT 
nO }0 T=t,NMAT : 
READCIN, 10°72) YMCE),PROT),GECI),DENSCI) . 
c YM=YOUNG#S MND./1u##h, PR=POISSON RATIO, GE=MUD. UF RIGIDITY 
C DENS=PENSITY : 
TECCECE EQN.) CECE VEYMCT 702.41. +PRCI))) 
YM(T)=YMC1)*#' FA 
© GECE)=GE(1)#1.Fo ; 
WRITECOUT, vic) CEOYMCT), PRCT) GECL) DENS(T), 1=L-NMAT) 
NO -bo PEt NMAT aa 
Poh HeENSCLI=SDENS(T)/0657.1/4412.) 
c INPUT JOINT COORDINATES 
REANCIN, 105.5) CXOMD> YOM), MAL NUTS) 
WRETEC OUTS > 04S) 
WRITECOUT.* 14) CHO XCMD> YOM) M=LONJTS) 
INPUT JOINT RESTRAINT CODE 
W=FREE 
1 =A LAMPED 
00 12 T=1,NJTS 
NRICE) =U 
NR? (120 
NRAC1)=0 
Nic) )=n 
N2(E)=8 
12 Nod] )=6 
TE CNR.EO.0) GO TO Ai 
HRITECOUT,’ HO) 
NY it T=t,NR 
READCIN, 19°51) STMT >M2, MS ; 


zm 


‘’ 


~ 
set 
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prec a aa er it 


11 
fea 


cae Eker Bea Bie Bi 


94 


qa°7 


Pos 


NRO CIE SM2 

NRSC ITIAMS 

WRITECOUT, u7) SToMILM2,43 

CUNT ENUF 

CONTINUE 

GENERATE COORDINATE NUMBERS FOR EACH. DEGREF OF FREEDOM, W IF 
CLAMPER, NOKMAL DISPLACEMENTS ARE NUNBEREY FIRST 7 
N}, N2, NS CONTAIN COORD. NUMBERS FUR EACH JOIN] i 
NREDU = NO. OF NORMAL NISPLACEMENTS 


CALL COORDNCNRI »NRP -NRSSNL NZONS ONIIS,NREDU, NDS ). 
HWRITEC OUT, un) 

DO 5 Jsl,NITS 

Tt CNRidd).6&Y.1) GO TO 8 

WRITECOUT, 2021) T,NI4E} 

CONT ENUE 

TNHPUT LUMPED MASSES s 
TEC NLUMP.EQ.f) GO TO 250 

READCIN, Tun) (CJMASSCT).RSMASSCE)). F=1>-NLUMP ) 
WREIECOUT, 27) 

Au -St T=L.NLUMP 

WRITECOUT, 925) SJHASSC IT)» RSMASS(T) 
RSMASSCLIARSMASSC19/ (5921/4812. ) 

CONTINUE 


NDF = NO. OF DEGREES GF FREEROM INCLUDING RUTATIONS | 


> CONTINUE 


NSSTFENDFaeCNDF At) 2 

NnQ 13 LT=l,NSSTF 

SST# C1) =0. 

TR CNBE.EQ.") GO TO cay 

REGIN TO GENERATE BEAM STIFFNESS TERMS 

WRITE COUT, 88) 

NO 14 NMFL,NBE 

INPUT BEAM FLEMENT PROPERTIES 

READCIN, TH 4) ARLXETSYISMAT.STNR, JTF R 

AK=A4REA OF BFAM CROSS SECTION, XI=SAREA MOMFNT OF INERTIA? 
YJ=EFFECTIVtE TORSIONAL MOMENT OF JNERTIA, MAT=MATERIAL CUDE 
JINR, UTFR=JUINT NUMBERS AT ENDS -. ’ 
WKRITECOUT,>U9NS) NMSARS XS. Y¥ So MAT. JINR, JIFR 

SkET UP CODF NUMBERS 

NOSE CLISNICUTNR) 

NUSUC(7 SNA CUTNR) 

NOSE (5 SNOCUTNR) 

NOSC(4)=N1(SUTFR) 

NUSG (>) SN2C ITER) 

NOS. (6 ENOCUTFR) 

TECMKEY.EQ.t) GU TO 25 

SINRE INFO. FOR LATER USE 

WRITECIDISC) ARsX1,YJ,MAT. JINR, JTHR, CONSE v1>0) 
CONTINUE 

XT=X¢CJSTNR) 

X¥2=XCSTFR) 

YIsYC TNR) 
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NRICJST)=AL | 


3 
pon we nae ements we aman 
«Ree Ue thas att 2, 


i ere “ 
Weis eee cate ee aes 
fi 


Ye=VYCITFR) 
FLNTH=SQRT((X2-X1 )@a2e (Y2-Y¥1 dae) 
GALL TRANS(X1,X2,Y1,¥2,FLNFH. DCS) 
o FOYM{ MAT) 
: G=CE (MAT) 
CALL BEAMKCELANTH.EG>X1,Y¥S+STM-DES) 
s Hi 96 Kets 
TF CNOSC(K).EQ.") GO To 15 
P=NUSCOK) 
NO 16 Nei, 
iF (MOSC(N).EQ.*) GO TO 16 
ASNUSE CN) 
Ts (J--LT.3)} GO 10 56 
MM= (2 e5¢(E-2 eC 2eNDF-1))/2 : a3 
SSTE CMM) SSSTF (Nr) tSTM(KSN) ! 
th CONT ENUE 
i& CONTERUE 
14 CONTINUE 
240 CONTINUE 
TE CNPE.EQO.')} GO TO on 
c REGIN TO GENFRATE TRIANGULAK PLATE STIFFNESS TERMS 
WRITECOUT,' nin) 
nU 17 NMS). NPE 
C INPUT TRIANGULAR PLATE ELEMENT PRUPERTIES 
RE AI iN,1t »5) PTH, MAT, .J71, JS? >, JTS. NDX 
PTH=PLATE THICKNESS,» MATS=MATERIAL CUDE, 
I11,3T2, ST‘ JOINT NUMBERS Al CORNERS, ANGLE 47 J35 MUST NUT BE 
8» LEGREES 
Nx WY, Dt DXY.HFTA - FLEXURAL RIGTUITY TERMS AND ANGLE OF MASERTAL 
PRINCIPAL AXFS wW/u FREANGLE LOCAL AXES 
TFCNDX.EQ.') REANDCEIN. 10027) DX,0Y,015DXY+BETA 
fFCNDX-EQ.!') GO TO 18 
RETA=u- 
DYSCYNCMAT I OPTHOS IDS CLP HC Tee PRIMAT I ##27)) 
nyY=)X 
Ni=PRCMAT) DY 
AXV=(C1.°PROMAT) 72. )9#NX 
18 RETASBETA/ (7.279598 
WRETECOUT, Ou11) NMOPTH MAE, JTL JP JTS ORXDY. UIs XY,RETA 
C; SET UP CODE NUMBERS 
NOSE CLI SNE CITI) 
NOSG C2 ENE CSTI) 
NOSLOS) ENS (UTI) 
NOSL(4)=NLC UT?) 
NUSL (9D =ENe( ITO) 
NOSC (6) SNSC UT?) 
NOSCC/)=N81 (CUTS) 
NOSC CH) SWZ CIT») 
NOSE CY) SNS CUTS) 
IF (MKEY.EQ@-1) GO TO 204 
C SIMKE INFO. FOR LATER USE 
WRITECIDEISC®) PTH. MAT JT1,3T2>J37o CNUSC( I), J51,9) 


ah 


NQMANAAIO 


A234 


7? 


CONT ENUE 

Bxi=X(JT1) 
RXP=K( ITZ) 
BYT=VCJTI) 
RY?P=Y(JT2) 


Yeo SORT (BX e-AXE Deas (BYZ-BYT ) #82) 

CALL TRANS(BX1.BX2+BYLeBY2,17,NCS) 

XS=NCS(2 8CXCITIV-“BX1 P-VCS CELIAC YC STS I-BY1) 
YS=VCSCIIOCXCITSIABX1 DES CZI ACV CITI -BYTD 

CALL PLATEKCY?2.X3,¥S.BXsD7- DI, OXY, RETA. DOS, PLIK) 


no 19 K=i,? 


1FCNOSC(K).EN.@) CGO TO 19 


T=NCSCCK) 
nO 70 N=1,% 


TE CNOSC(N).EQ.9) GO TU 20 


J=NUSCCND 


IF CSeLT.D) GO TO ot 
MM= (Pes Ci-1beC2enDFe L))/¢ 
SSTF (MM )=SSTF CHM) © PL IK (K,N) 


CONT ENUE 
CONTENUE 
CONTINUE 
CONTINUE 


STORE FOR RENUCTION 

NO <<} Pst, NF 
NSe(PadeCh-1) eC ZenBF ef) )/4 
NE=(PaeNDF eC hf 8 ( CONDE 1E)D/2 
WRITECHBISC) (SSTF CJ), JENS. NE) 
REWIND IDISC 

NU 22? [si,NSSTF 


SMC 1)=0. 


TF (MKEY.EQ.1) GO TO Ao 

TECNBE.EQ.") GHA TO £01 

GENERATE HEAM MASS MATRICES 

NO cd NM=1.NBE 

REANCIDISC) AR, XEsVJIs MAT» JINR. JTFR. (NOSUC IT), F556) 


Ku=XCJTNR) 
X7Z=EXC ITER) 
YI=YCJTNR) 
¥¢CHVCSTFRI 


FLNTHSSORT (CX 2-X1 O62 + CY2~-¥1 Dee) 

CALL TRANSCOX',X2,V1e YZ -FLNTH- DES) 
RHO=DENS(M4{) 

CALL HEAMMCEL NTH. RHO, ARS KE YJS>SMM, MUS) 


na «4 K=3,¢ 


ib (NOSC(K}.ER.' ) GO TU 24 


T=NOSC(K) 
nO 7-5 Net,” 


TF CNOSC(N).£60.9) GO TO 25 


J=NOSCCN) 


IF(UeLTE) GE TO 45 
MM=(2ase( 1-1) eC 2enDF-1))77 
SMC MM) =SMCMM)D+SMMCK SND 


aes Eve eines tee 


al 
3 


25 CONTENUE 
“4 CONTINUE 
24 CONTINUE 
2ut CONTINUE : 
iF (NPE.EG@.:') GO TO sft 
C GENERATE TRIANGULAR PLATE MASS MATRICES 
DU 76 NH=1.NPE 
READCIDISC) PTH MAT, JTt,JT?,J57T5,(NGSECE), 221.9) 
AXT=X¢( STI) 
BXP=X( IT?) 
BYt=Y¥(JST1) 
By72=¥(JST2) 
Y7=SORT (CBX. -HXL) BOC (CRY?-BY1 #82) 
GALL TRANS(BX1.BX¢-BYi1.BY2,Y7,fC0§) 
XV2pOSe? deCKC STS )-BXUP-DCS(CL OCC ITS) -BYT } i 
Yo=pOS 1 .eCX( STS) OXI DCS (2d aC VC STO) HBYI D 
PRHU=NENS( MAT) 
CALL PLATEM(Y?.X5,Y¥5,PRHO,PTH,DCS,PLTM) 
N00 «7 K=1,% 
Tt (NOSEC(K).£Q.0) GO TU 27 
T=NOSC CK) 
fu ¢® Net," 
It (NOSC(N).E0.") GO TO 2&8 
JSNOSCON) 
TRC sel TH) GO TO 78 
MM=(7edeC]-Te(2enDF-1))/2 
SMC RM)=SK(MM) + PLIMC KN) 
chk CONTENUE 
7 CUNETANUE 
cf CUNTENUE 
tut CONEENUE 
C STARE FOR RENUCTION 
245 CONTINUE 
NO -Sd FFL. NLUMP 
NN=JMASSCI1) 
IF ONT CNN) FU.) BO 10 758 
NNN=N1 (NN) 
NS= (7 © NNN¢ CNNA-~ 1) 807 ONDE ONAN) 7/2 
SMONS)=SMONS)+¢RSMASS CHA) 
ok CONT INUF : 
nu .¥9 Tet, Nur 
NS= (Zadeh) Ved eanDbe())/ 7 
NE=(CAONDE eC 1-1 00¢ ONDE 2) 72 
oO URTTECNDISO) (S405), J3=2NS.NE) 
NUMASS=NDF-NREDU 
CALL FIGENCA. VALU. TEMP.H,C. UUM’, F,TDUM4, EbISC, JDISC,KUISC,NUISt, 
MU TSC NDE NMGNE, NMODE,NREDU,NOMASS) 
RO 10 100 
ENN 


~ 
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FORTRAN DECK . i 
TCOORBN ASSIGNS A CUORD. NO. 10 FACH HEGREE OF FREEUOM AT EACH JOENE 
NRILNR2,NRS = ARRAYS CONTAINING RESTRAINT INFO. FUR EACH DEGREE 
NF FREEDOM AT FACH JOINT (FREESu., CLAMPEN=1) i 
NI.N?,N3 = COORD, WO. FOR FACH DEGREE OF FREF NOK (NORMAL 
DISPLACEMENTS ARE NUMBFRED FIRST) 
NJTS = NO. OF JOENTS 
NREDU = NO. OF NORMAL DISPLACEMENTS 
NDF = TOTAL NO. OF DEGREES OF FREENDUM CINSLUDING ROTAITUNS) 
SURROUTINE COORDNCNRIDNR2 NRO NL ON? ONS /NJTS,NREDUSMUF ) ; 
DIMEASTON NRECSU}.MR2(50) ARSC Ud RLCOD) NPC 0) NOCOH) [ 
NOS} 
NO 30 [=1,NUTS 
IFCNRLCIT).FO.13 6 TO 10 
NiCTD=NO 
NU=NO+1 
CONTINUE 
HRFELUZNO=1 
BO ch Y=i,NNTS 
TF CWR2(1).F0.7) 86 TO 71 
NW? CL) ENO 
NO=NO+1 
21 TECNRSOC1).FO.1) GO TO 7U 
NOCT END 
NO=NO+1 
20 CUNTENUE 
NDF=NO-J 
RETURN 
EAN 
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is FORTRARK DFUK 
CTRANS TRANS? ORMATION DIRECTION COSTNES 


c Xt. Y1 = COMRNS. OF POINT } : 
c X7.Y2 = COODRNS. OF POINT 2 ; : 
'y FL = DISTANUF BETWEEN POINTS { AND ¢ ; 
C NCS = DIRECTION COSINES OF VECTOR FROM POINT i 10 PUIKT 7 4 

SURROUTINE TRANSEX!>X2.¥1,%2sFL.0US) d 


NIMENSTION NUS (7) 
NUSELSCK2S-K1 FL 
ACSOP VECVAWYU FL 
RETURN 

FND 


tad RATER ot, 


f FORTRAN DECK ee 
CHEAMK PLANE GRID BEAM ELEMENT STIFFNESS MATRIX IN SYSIFM LOUKDS. 

c F) = BEAM LENGTH 

¢ & = YOUNG*S NONULUS ‘ 
c GR = MODULUS OF RIGELITY can 3 
c ¥1 = ARGA MUMENT UF INERTIA : 
c Yj = EFFECTIVE TORSIGNAL MOMENT Gt INERTIA 
C SIM = STIFFNESS MATRIX ay | 
c NES = DIRECTION COSINES 


SURROUTINE BEAKK(F LE, 62X1+ Yds STM, DUS) 
NTMENSTON S1M(6,6),0CS(2) 
ZiZbOXIT/FL 
72208 ISFL 
SIM(1,T)BL oe ZI CFLAFL) 
S1H(2.1) 26. #71 8DES(2/FL 
STM( 2.2984. #71 HDCS(Z ANCES?) +27 ANUS (TD eDES(L) 
SIMC3, 12) Soe eZ DESO) FL 
SIMS 2 DECK 4 A ZEE LZ DOMCSCL MDE CA} 
SIM(A,5)=4 U7 FE MDCSCL ADESC F274 DESC? EDESCA) 
S1M64,) =~STM(1e?) 
S164. 2PEeSEMC 7,F) 
q STH04,5)=-STMC 401) 
2 STMC4.4) ESTA bd) 
4 s1M(5.1)SSTMC2- 2b) 
SIM, 7IEZ- #21 0CS(7 ONES (2 I~Z2 4 DCSE DESC) 
“fl STM(S 3) E— C7 #7 L4Z2 ODESCL ODES (2) 
: STM(5,4)E<STM( 701) 
q S1M¢5,9)5STM(2-2) 
STH(6,1)=STM( 5-1? 
SIM(H,2 F5THOD6 5) 
4 STH( 6, 3)22- #71 #0CSEI NES (1 -L78DUS(7 DHBES Ce) 
f STM(6,4)5-57M( 69) 
f SIM(6,5)=STM(5-2) 
S1M(16,6)=SIM(0S+5) 
no 1A ps2. 
N=]-1 
no it Jel.N 
yt STMCd, PD =STMCL J) 
RETURN 
END 
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$ FORTRAN UFCK 
CHEAMM PLANE GRID BEAM FLEMFNT MASS MATRIX FN SYSTEM COOKBS. 
C FL = BEAM LENGTH 
c RHO = DENSTTY 
C A = CROSS SECTIONAL AREA 
C XI = AREA MOMENT OF INERTIA 
tf XJ = EFFECTIVE TORSIONAL MOMENT OF INERTIA 
c Sm = 4aSS MATRIY 
t AES = DIRECTION CUSTNES 
SURROUTINE KFAMMCFL > RHODA,X3-7%75,SMM>DUS) 
~ DEMENSTON SMH(4,6),BES(2) 
ZiZRHOSASEL 
7eztl ee? 
7I2KI7A 
NPs Ptah 3. 13564 (6.07597 (9.422755 
CCr Zeca estth ss’ td, 473/010. 8FL)) 
AAzZVE(CZA/ Uo eter aZs/ ihe) 
TH=Z18Y¥U/ CS. A) ’ 
RR=Z1 (Dos ‘to lo.#Z0)/ (5,822. ) 
QUSZIEC VSL SIAN. HZ5/(100.0F RD) 
SSr-Z1e(Z2 pt. eZo7a,) 
PP=Z1TSYIS(%. 8A?) 
SMécd, 1 =00 
SMM(?, ft FCC ENS (7) 
SMM(2,27 DEAAMNOS(C? OBES (29 +T1 #008 (1) # DCS (1) 
SMMt4.7 2 F-CCe NCS) 
SMM(5,27)2(-AASTT) aDOS(1 DeUCS (2) 
SMM(.3,5)2AA@NES( 1 OVES(1> + TIDES (7) #DCS(2) 
SMM(4,17)=RR 
SMM(4.?)=ON#NCS(?) 
SME14, 4)=-NU8PTSC1) 
SMM(4,4)=SMMC().L) 
SMN(5,12=-SMM( 4,7) 
SMM(5,7)=SS5eNES C2 PVCS (Ce d+ PPHHES C1 peDUS(t > 
SMM(5,3)=(-SSe+PPIEDES (1 eNOS (-) 
SMM(4%,4)=-SMHM(",1) 
SMM(%,5)=SMM(27.2) 
SMH(46,1)=-SMM(4, 3) 
SMM(6,2?)=SMM(. 5) 
SMM(6,3)=SS#HCOS( i ODCS(1) + PPeLDCS(, >*#BCS(2) 
SMM(6,4)5-SMM(551) 
SMM(16,59)=SMN( 5.27) 
SMM(6,5)25NN( 5.5) 
hu if [Ts2.% 
N=f-1 
nO it J=1,N 
1h SMMCJS,T)=SMM(]. Jd) 
RETURN 
FANN 
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THIS SUBROUTINE: DETERMINES THE STIFFNESS MAIRIX OF A 

TRIANGLE PLATE ELEMENT IN SYSTEM COORDS. j 
Y7.x3,YS = COORDS. OF PLATE CORNERS IN LOCAL COURDINATES : 
OX, UY,DI+UXY.BETA = FLEXURAL RIGIUITY TERMS AND ANGLE O¢ MATEREAL 4 
PRINCIPAL AXFS W/O TRIANGLE LOCAL AXES 

HCS = DIRECTION COSINES 

PLTK = STIFFNESS MATRIX 

SURROUTINE PLATEKCYZ,X30YS/ DK, DY, 01 -DXY+RETA,DCS»PI IK) 

DIMENSTON PLTK(9,9)26(959) 2 CTNV( 99) 2 P(9,9) oR(90Y) 

DIMENSION T(4,9)/STIFF(9.9),DCS(2) 
EQUIVALENCE (P(1,.: )+STIFFCL22))> CRUDE ECL dD? 

CALL CMATCYZ. XS V 500) 

CALL MINV(C,CINV,Y) 

CALL DENMATCY2.X3,Yo,DX,DY.D1,DXY-RETA>P) 

CALL MATMPY(P,CINV.R,9) 
n0 10 [22, 

N=I-1 

NU 10 J=1,N 
ZZ1=CINVOIL- J) 
ZZ?=CINV(J-T) 
CINVC EL, J)=22? 

CENVGJS, T2221 

CONTINUE 

CALL MATMPY(CINV.R STIFF.9) 
nO 4hy Je1.4 

nO 40u Jel.y 

TCT, Sate 

T(1T,1}2). 

T(4,4)21. 

T(7,7)=1- 

T(?,?)=2DCES(z) 
T(%,3)=NCS(2) 
TO4,5)=DCS(2? 
T(6,6)20CS(2) 
TCR, 8) =DES(2) 
T69,92NCSle) 
T(?,3)2-DES(1) 
T(5,6)2-DOS(1) 
T(R,9)=-DCS(1) 
TO%,?)2DOS(1) 
T(46,5)=0C0S¢1) 
1(9,8)2DCS(1) 

CALL MATMPYCSTIFF,T.C.9) 
T(?,3)2DCS(1) 
T(h,6)2DC5(1) 
7(8,9)=DCS$0)) 
143,2)2-DCS(1) 
1(6,5)2-DCS(1) 
1(9,8)2-DCS(1) 

CALL MATMPYCT,C>PLTK.Y) 
RETURN 
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FORTRAN DECK 
CHAT ss 
THIS SUBROUTINE FORMS THE C MATRIX RELATING THF CURRER 
HISPLACENENTS TO THE POLYNOMIAL DEFLECTION COEFFICIENTS 
FOR THE TRIANGULAR PLATE ELEMENT 
Ye. X3,Y5 = COORDS. OF PLATE CURWERS IN LOCAL COORDINATES 
mos C MATRIX 
SURROUTIME CMATCY7>X3,Y550) . 
DIMENSION €(9,9) 
NGO i [21,% 
nu 10 Jet.” 
18 CCT,JI=0- 
Cli,F RL. 
CP +-3aL.) 
C(3,2)2-1. 
C4, zh 6 
C(4,3)=Y¥2 
C14,6)2Y2 en? 
AA, 9)a¥Paas 
C(S,3)216 
C15,6) zee #¥Z 
C05,9 )=5e8 VOR? 
C(A,? pst. 
06625) 2-V¥¢ 
C(6,8)=-¥eue? 
C(7,1)=16 
6(7,2)=X3 
0(7,3)=Y3 
CO7 4) EXSae? 
0(7,5)2XdeYa 
O07 .6)5¥ See 
COl.7=XS HF 
GO7,B)=KIaVSeareVsaxseaes 
(7,9) )2V¥Saed 
OCR, 3)=Le 
CCR,5 2X3! 
CCP, 6)=208VO 
COR,ADES ce AXSTV SX See? 
CRISS HV Gea” 
C(9,? pz-t, 
6(9,4)=-2,4X% 
(9,5 )2-¥5 
C(9,7)=-3.#X3aa2 
C(O, Bp =HLY seevee. aXveys) 
RETURN 
END 


——¥ 
3 


HRI A 


$3 


nn nemeintiinenatiianan tiie to saieceraddraassroseseancnatenns caiman anemmencihemtamnemmnnaeamiadneindiaimennh mamas eemen aimee aon 
° 


© eng ew 


t FORTRAN DECK 
CKINV MATRIX INVERSION SUBROUTINE 
e A MATRIX 10 BE INVERTED 
¢ U INVERTED MATRIX 
c NM = ORDER OF MATRIX {.LE.9) 
SURROUTINE MINVQAsU2NM) 
NEMENSTON A(9,9),U09.9) 
no yfut T=1,NM 
no y0ul J=!,NM 
UC]. J seb 
TF (CE.EQ. J) UCI,J)22.0 
gnnt CONTINUE 
EPS=0.NnQutoN) 
ne vyA1ts F=', AM 
K=1 
TF (T-NM) 9U71-.9007-90271 
gngt TF CACHE» T)-€PS) On05, 9NV6,90"/ 
yhusS Th Coad LP-EPS) 906,906. Ynu7 
gnns K=K+) 
pO ¥t23 J=!,NM 
UCT» J)Sudd. JdeUCK dD 
gnes ACL JZACE JDeACKS SD 
GO 10 921 
gniiz Niveath,T) 
NO vud J=),NM 
Ute J>euch Jd7DLV 
yang act, JZaCl J/DL¥ 
00 ¥O0L5S MM=1,NM 
MEL T=ACMM, 1) 
If (CABS(DELT)-FPS) Y01L5,901D VELA 
ynia TE CHM-7) 694,901 5,90TU 
vita fa YOlt Je) oe NM 
CMM, JPSUCMM. JP-UCE JDODELT 
gor AC MMs JDEACMM. JDACT, J) @DELT 
gai CONTINUE 
00 vss Tet. NM 
no ¥y053 Jz) NM 
ghod ACE s)eUCT dS) 
RETURN 
FENN 
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¢ FORTRAN DECK 

COINMAT 

t THIS SUBROUTINE DETERMINES THE DOURLE INTEGRAL MATRIX FOR. 

¢ THE K EQUATION FOR THE TREANGULAR PLATE ELEMENT 

C ¥2.x3,745 = CHORDS. OF PLATE CORNERS IN LOCAL COORDINATES 

t DX,UY,D1l,DXY.BFTA = FLEXURAL RIGIUITY TERMS AND ANGLE OF PALENUAE 
¢ PRINCIPAL AXES 4/0 TRIANGLE LOCAL AXES 

C P = DOUBLE INTFGRAL MATRIX 


SURPOUTINGE UINMAT(CY2,X3,YS,DX.DY,U1+DXY>BETA>P) 
DIMENSION P(9,9).D00,5) 

00 3 Ji, § 

nga yh J=i,9 

P(1,J)=0- 

CALL HHATCDX.BY,01-DXY,BETA.D? 

AL=DBLINTCYZ. X32 Yoo od) j 
AZ=BBLEIENTCYZ.X3,VSe1,0) : 
ASZ=PBLINTCY?. XA, Yo,e00) 

A4=DBLINTCY?. X3,VSsto Ll) 

AD=UBLINT CY 2. X36 ¥o0 be?) 

AOZVBLENTCY?. X4eVoodo Ll) 

P(4,4)=4-4041,1 0A) 

P(4,5)2 4,.4DE1,5)0AL 

P(4,6)545.4N(1,7) #Al 

P(4,7)=12.#D(1,1L) #Az 

P(4,8)=4oR( NCTE LIHASHDC1, ZdPAS+ eZ ADC TT) aC AZ+A4)) 
P(4.9)212.40(1-2)0A4 

PC, 5)=4o#N( , 5) 01. 

P(%,6)> 4.2064.-2) 847 

P(&,7)= 12.8451 eA? 

PCH, BY 468(NC4,1 )CASHH(S. 2 2HAP HP ODI Sr» SIOCAZ+A4) ) 
P(5.9)= 12.8N1(.5,7) 4A4 

P(6,6)=4—411(7,? Daal 

P(A,7)=12 64h (21) ah 

PCA BRM H(BC QL PAVED( Zr 2d HALES eo ANC As S)BCALIAY)) 
P(659)=126#D(7,2)8A4 

P(7,7 SSG CL LIMA’ 

PC7T/BYZL2Ze (AC 1 1 dHAGEHC1 2 dFASEA LUCA SD) ECASEAO)) 
P(/,9)=36.20(1,2)8A6 

PCR, B)=4e#( DCT LIBPADFEN (1,2) #ADHzS ODI, 3S) eCAGtA*)) 


1 +4oF( D2 LIBAOFD (2,2) HAS+?2 eer S)aCAdIAS)) 
1 tHe t(D Se LD eAOt( Ss 2dHASHEZ.ONIS, SO CASHAD)) 
1 Bet (NCSL VOAD*E (Ss 2d HADFHZ ANI S, STS) B(CAOH AD) 


P(R,D)=R7AoH(DC 12) OASFHD(2,2 MAKE OD S22) BC ABAD)? 
P(9,9)=56.#D( 2.2) 8A5 


00 70 [=1.,° 

Ne}el 

No 7-0 Jen,” 
P(J,T)=PCE. J) 
RETURN ‘ 
END 
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¢ FORTRAN DECK 

CHATMPY 

Cc MULTIPLIES MATRICES A AND B TO GET C, ALL OF ORDER NON 
SUBROUTINE MATMPY(&,8,0,WN) ‘ 
DEMENSION A(9,9),809,9),0(9,9) : 
ne 16 Jat,n i 


CCL, J ate 
NO 10 Ks1,N 
18 CCT, SACL. JAC KIeB(KSS) 
RETURN 
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t FORTRAN DECK 

CASLIRT ; 

c THIS SUBROUTINE EVALUATES THE DOUBLE INTEGRALS APPEARING IN THE 
c EQUATIONS TUR K AND M FOR THE TRIANGULAR PLATE ELEMENT 

c Y?.X3,Y3 = COORDS. OF PLATE CORNERS IN LOCAL COURDINATES 

c M,N = POWER OF X AND Y RESPECTIVELY, PRZEMIENIECKI, PAGE vtld 


FUNCTION DRLINTCY2+¥3,Y3>MeN) 
DIKENSETON AX(2)-816-),P2073.P2(7),P9¢7) 
FOUIVALENCECRILGCL).PSCL)) 
TF (Mm1) 4u,-41,42 
40 P)(i)=1.6 
NI=b 
60 10 45 
41 PillIs-Led 
Pi(e)=1.0 
Niz=j 
80 TO 43 
4? CONTINUE 
Alb(Qadz-1eU 
Al(c321.0 
Qy(l pst 
Bl (2 )s1.0 
M11 
44=K~-i 
nO 30 J=2,MH 
CALL PLYMP(A1,1,81.M1,P1,N1) 
NNT=N1 41 
HO 30 P=1,NNt 
Ai td=PI(7) 
MI=N1 
19 CONTINUE 
43 CONTINUE 
EF CN-1) 50,51,52 
5N PZ(2)=1.9 
N2=U 
GO TO 553 
61 P21) a-¥S+YV? 
Pee )AYS 
N?=] 
60 10 55 
5? CONTINUE 
AIC) =-Vd*¥e 
Ai(c)=¥3 
RVC) V=E-¥S+V? 
Ql cdaY3 
412) 
NN=N-] 
NO 70 J=1,NN 
CALL PLYMP(A1,1,.B1241,P2.,N2) 
NNP=N241 
nO ¢ J[21,NN? 
BLCtd=P2C01) 


58 


Mi=N7 
70 CONTINUE 
23 CONTINUE 
CALL PLYNP(P1,N1L.P7,H2,P3,NS5) 
WNS=N542 
Sol=t. 
nO 30 F=1,NN% 
SOLESOL+ (CX Soe (+t eV 2a LT. SELUATOMOMNOZ) De PECL ACES SFLGATINS+2~1)) 
aM CONTENUE . 
rot INT=SOL 
RETURN 
En 
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q FORTRAN BECK 


CPL YHMP L'7e§ 
& POLYNOKHEAL MULTIPLY 
SURKOUTEME PLYMP(A,Loa-6- C,H) 
C Cl-3ue Quim-hAg 
mt HEMENSION ACLEoB(I).00L) 
| ¥ N>Le Wf 
{128 4+ |: 
OG HT = «obs 
u a i ee) 
LF =t¢#: 
4? = + 
bo} t 2 ,t? 
On 1 J = ',M? 
K = [6d 
1 CiKei) = CCRH19 s aC bseB( dS) 
un ? K = oth 
it = K 
1FCCCKI) 4.7,3 
CONTINUE 
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J 
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FORTRAN DICK 


PLATEN 


4an 


THIS SUBROUTINE NETERMINES THE MASS MATHIX UF A 
TRIANGLE PLATE ELERENT IN SYSTEM VOURDS. 
¥e.xX3,%09 = COORDS. OF PLATE CORNERS IN LOCAL COURDINATES 
PRHOG = DENSITY : 

Pik PLATE THICKNESS 

nes DERECTION COSINES 

PLTM = MASS MATRIX 

SURROUTINE PLATEMCYZ. X45 %9,PRHO,PTH,DUS,PLTRM) 
DIMENSION PLTM(9,9),0(9,9) CINV(9,9),P(9,9),R(°9,¥) 
DIMENSION 7(9,°),FMASS(9,9),00S(?) 
FOUIVALENC! (P€1,1)-FMASSC L512), CRCt,9)-TC1,45)) 
CALL CMATCY? .X3475,€0) 

CAIL MINVCC,CINV,Y) 

CALL DINMTHCY2Z.X3,Yo.,PRHO,PTH,P) 

CALL MATMPYCP,CIANVOReY) 

HO 3: Js2,%, 

N=1-7 

NU IH Jet. 

771=CINVOCE.J) 

7TL°O=CINVOJS- ED 

CUNVC EL, J)e22? 

CINVOJ,F=7Z! 

COMTERUE . 

CAEL MATMPYC(CCINV,ROFMASS.9) 

Ng «flu Jebey 

M- aft Jtr-y 

TCL, Jdsu. 

TCT, T)=be 

T(4,4)216 

T(7,7)=) 6 

T(?,? =VES CY) 

TCS, 8) MOS 2) 

TC, 5) SNCS( 7) 

Tlo,4)20CS(2) 

Fi, 8)=VCS(7) 

T(O,9) SNCS. ) 

T(7.3)2-DCS(1) 

1(5,6)=-DCS(1) 

TCK,9}=-DES(1 

T(%,2) DESC) 

T(4,5)2nC0S(1) 

T(9,8)=NCS017 

CALL MATMPY(FMASS,7,C0,9) 

TO?, SENCS(1) 

TO5,6)20CS(1) 

TCR, 9 EVES (1) 

T3627 )=-BUSC' ) 

T(6,5)=-DUS( 7 

T(9,8)2-DCS(1) 

CALL MATMPYCT.LPLIM.Y) 

RETURN 

END 
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+s fOKTRAN DECK 

COTAMTM ‘ 

c THIS SUBRQUTINE DETERMINES THE DOURLE INTEGRAL MAIRIX FUR 
c TRE TRIANGULAR PLATE M MATREX - PRZEMIENTECKI, PAGE 304 

ys Y¥eoX3,Y%o = COORDS. OF PLATE CORNEKS IN LOCAL COORDINATES 
t PRHO = DENSITY 

c PTH = PLATF THICKNESS 

c P = DOUBLE INTEGRAL MATRIX 


SURROUTINE DINMTMC Ye. X3,Y5S,PRHO,PTHOP) 
DEIMENSTON P(9,9) 
PCT,TISDBLINTC V2. X3,V¥5.4,0) 
PC?, 1) =DBLINTCY2,X3,¥5.1.0) 
PC7,?2 ZSDBLINTCY2,X3,YS.250) 
PCS, PISDBLINT CYS]. X3,V500>1) 
P(4,2 =NBLINTCY2. X5,Vaed> 1) 
PCS, SPEPBLEINT CY» X3,V¥5,U,7) 
P(4,1)=P (2.7) 
P(4,27 }=NBLINT(Y2,X3,¥3-3,0) 
P(A, SIZNBLINTCYZ X35 ¥50251) 
P(4,4)2DBL INTC V2. X3,V50450) 
P(S,1)=2P05-¢) 
P(5,2)=P04.5) 
P(4, SSDBLINTCY2.X3,¥5.15¢) 
PCS, 4) =NBLINTCV2>X3,V5,551L) 
. PCS, 5) SDBLINTCY2.X3,%92252) 
PCA,TIEP CO.) 
P(6,7)=P 09-0) 
P(A, S)ZADBLEINTCY2,X3,Y9-U,5) 
P(4,4)=P(d.5) 
P(6.5)=DBLINTCY2-X3,Y5-259) 
PCA H)SDBLINTCY2sX3,Yo+0,4) 
P(7,1)2P04.7) 
P(7,2)=P(4.4) 
PCT, )=P(b.4) 
P(7,4)=DBLINT(VY2-X3,V5.5,0) 
P(7,9)2NBLINTCY2,X3+¥9.451) 
P(7,6)=DBLINTCV2-X3,V¥3-5,2) 
PC7.7)=NBLINTCY2.X4,V3.0,>) 
PH, 1)=P (0.594004, 4) 
P( 8,2? )=P(0.49)9+P(5,4) 
P(H,3)=P(0.57¢P (5,5) 
P(R,4)=P¢0/,.60)4P(7,5) 
PCR, 5VEDBLINT(CY2.X43,15.2553)%PC7,6) 
PCH,O)SOBLINTCYZ.X3,Vi- LF tDWL INT CV 2,X5,VSe-, 1) 
P(R,7)SDBLINTCY2. X3,V5,4-23¢DBL INTC VY? XS,VG5,') 
PCR,A=DBLINTCY2. X35, V5 e254)+2BLINICV2,X3,¥S,9,7) 
1 +P. H#DBLINTCV 2. X35 V¥5555 9? 
P(9.1)=P(0.9) 
P(9,7)=P(0.5) 
P(9.45)=P (6.6) 
P(2,4)=PALINTC V2). X4-Ya.25 5) 
P(9,5)=DBLINT(V2.X3,¥9,1,4) 
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P(9O,6)=NBLINTCY2. X3,V4:0.7) ‘ 
P(9,7)=DBLINT(¥2,435Y5,55-5) : ; 

P(OY, BYZUBLINT( V2, X3,Vi21,9)¢DBLINT(Y2,X9,1%55754) 
P(Y,9 )=DBLINT(Y2,X3,Y38,U54) 

nO 30 Jsi,” 

NO 40 Jzi,i 

P( J}, J)=PC1. J) *PRHUSPTH 4 
NU 20 T=2,%' E 
N=]-1 3 
00 70 J=1,N 

PCG, TISPCL. J) 

CONTINUE 

RETURN 

FND 
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FORTRAN DECK 


REDUCES STIFFRESS MATRIX AND INVERTS IT, REDUCES MASS MATRIX 
DETERMINES EIGENVALUES AND EIGENVECTORS 
THE ARGUMENTS ARE= 
A - VECTOR OF LENGTH NROFYCNRDE +1) 7/2 
VALU - VECTOR NF LENGTH WEICG 
TEMP.8,C-DUNt. - VECTORS OF LENGTH NRDF OR NBASS (SMALLER) 
E - MATRIX OF DIMENSION (CNRDF,3) 
IDUM4 - VECTOR OF LENGTH NRDF OR WHASS (SMALLER) 
ITAPE, JTAPE, NTAPE, MTAPE, - THESE ARE VARIOUS TAPES 
WRDF - NUMBER OF DEGREES OF FREEDOM OF THE SYSTEM 
NEIG - NUMBER OF EIGENVALUES DESIRED 
NVEC - NUMBER OF EIGENVECTORS DESIRED 
NMASS=NO. OF NORMAL DISPLACEMENTS 
NOMASS=NO. OF ROTATIONAL DEGREES OF FREERON 
STIFF 1S ON MTAPE IN COMPACT FORK 
MASS IS ON NTAPE IN COMPACT FORK 
SUBRGUTINE EIGENCA, VALU, TEMP,.8,C,DUM3.E> 1 DUM4, 1 TAPE» JTAPE,KTAPF, 
INTAPE, MTAPE, NRDF.NEIG,NVEC.NMASS, NOMASS) 


BIMENSTON DUMS(NRDE). 1NUM4(1),A(1), oe ee C(1),EC(NRDF,3), 
LTEMP (1): 


DIMENSTON ELOW(50), THIGH(59) 
INTEGER OUT 
MATA Q1/6HFLEXIB/-Q2/6HILITY /.Q3/6HMATRIX/ 
DATA O4/6HWEIGHT/.O5/6H MATRE/.Q6/6HX / 
NO 56 TI=1L.NMASS 
TLOWCLT) 23 
56 LHIGHC(IL)=NHASS 
NUT=6 
REWIND MTAPE 
REWIND NTAPE 
NTENP=NMASS 
CALI. DIVIDCNHASS, NOMASS,MTAPE. JTAPE> ITAPE. A,B) 
CALL ZROMAK( A,B,C. DUM3.NMASS.NOMASS, ITAPE, JIAPE,MTAPEKTAPE) 
CALL DIVIDENMASS, NOMASS, NTAPE, JTAPE. ITAPE, AB) 


CALL ZROMAM( A,B,C. DUNS. NMHASS.NOMASS, FTAPE. JTAPE,NFAPE, KTAPF) 
345 CONTINUE 


REWIND MTAPE 

REWIND NTAPE 

NREDUZNMASS 

NRMX=NREDU® (NREDU41 7/2 
READCHTAPE) (ACT), F21,NRMX) 


EISEN 


AHA VGA GaAaIaqarinrrrc 3 cr 


WRETE (OUT, 5590) 
5509 FORMAT(///65HR E DUC E D UPPER TREANGULAR 
1S TIFF WES S ATR I X) 


0 5507 [=1,NREDU 
NS=(2e]74(1-1)e(2¢NREDU-1))/2 
NE=(7eNREDUt (1-1)9(2@NREDU-1))/2 
WRITE(OUT,55"2) 1,¢A(5),SENS.NE) 

559? FORMAT(C/SHROW, 14.4(0/9E1425)) 

5591 CONTINUE 
CALL SYMINVCA,NREDU) 
WRITE (OUT, 5593) 

550% FORMAT(///69HR E DUCED UPPER 

irFLeXxTBTrEtEITY MATRIX) 


PUNCH 5662, CCTLOW(CK). FHIGH(K)),K=1>NREDU) 
560? FORMAT (1814) 


nO 5504 1=1,NREDU 
NS=(2e Te 1-13e(2eNREDU-]) 372 
NE=(2eNREDUt(T-L)¢¢2@NREDU-]))/7¢ 
5504 WRITE(OUT. 5597) 1,(ACJ),JS2N6.NF) 


TRIELANGULAR 


104 


611 


559 
5508 


551" 


6nin 
5507 


560 


570 
58 


691? 


5516 


5506 


5513 
5512 


5514 


5511 


355 


wWF TIGHT MA 


PUNCH 6011, 61,07,-02 

FORMAT (3A6) 

NO 5597 1=t,NREDU 

l= }et 

IFCST.EQ.9) GO TO 5508 

NO 5509 J=1,11 

NUs(2efe(do119(201-5))/20¢J=1 0 (NREDUOT) 

AC J)ZACNU) 

CONTINUE 

NS=(Pele( lo] )e(2eNREDU-1))/2 

NE=(2eNREDUS (J -1) eC 29NREDU~13)/2 

Jet ‘ 

NO 5510 JJ=ZNS,NE 

ACSA SI) 

Ja Jet 

PUNCH 6010, (R( J), J21,NREDU) 

FORMAT(1P6E12.5) 

CONTINUE 

OPTION TO EXPAND REDUCED FLEXIBILITY MATRIX TO FULL MATRIX BY 

INSERTING 1 OR 2 ZERO ROWS AND COLUMNS REPRESENTING ATTACH POINTS. 

CONF , NCOD = 1 OPTION EXECUTED , NCOD = 0 OPTION NOT EXECUTFD 

REAN(5,560) NCOD 

FORMAT (16) 

IFCNCOD) 580.580.5780 

CALL FULFL (A-NREDU) 

READ(NTAPE) CACT). F21,NRMX) 

DO 60412 I=1,NRMX 

ACT)=ACE)032.174812. 

WRITE (OUT, 5595) ; 

FORMAT(///79HR E DUC ED UPPER TREANGULAR 
TRI X) 

DO 5506 T=1,NREDU 

NS=(Pat4(To1)e(2eNREDU-1))/2 

NE=(2eNREDU?( 1-1) 84 7@NREDU-1))/2 

WRITECOUT,5592) 1,(AC J), JENS.NE) 

PUNCH 6011, @4,05,06 

NO 5511 T=1,NREDU 

list 

IFCTL.EQ.0) RO TO 5512 

BO 5613 J=1,71 

NUS( Pets (So1 (281-5) )/246(J91 )*(NREDUMT) 

ACH SACNU) 

CONTINUE 

NS=(PeleC1-1)9(2eNREDUH-1))/2 

NE= (2 aNREDU+ (1-1) 8(2@NREDU=1))/2 

Jz! 

NO 5514 JJENS,NE 

BC J) ZAC IS) 

J=J+1 

PUNCH 6010, (Q(5).J21,NREDU) 

CONTINUE é 

IF (NEIG.EQ.0) RETURN F 

NMAX=NTEMP®CNTEMP+1 )/2 

CONTINUE 

REAL IN THE MASS MATRIX 

REWIND NTAPE 

READ(NTAPE) (ACT), 121, 0RAX) 

REWIND NTAPE 

CONTINUE 

CALL EIGMATCNTEMP. A, VALU, TEMP,R,C.DUHY.E, (DUNG, NYSEESNTAPEDUTARE, 

TETAPE,NETG.NVECG) ; 105 


108 CONTINUE 

DO of f=1,NEIG 

DUMS CPI=SORT(VALUCID)/6.2831853 
60 CONTINUE 


WRITE (OUT.9099) ; 
WRETECOUT.90N5) (C1, DUMACT). 1=1, NETS) 
9909 FORMAT(1H1,43X.S3HHERE ARE THE NATURAL FREQUENCIES 471) 
9905 FORMAT(35X.,29HTHE NATURAL FREQUENCY WUMBER 15,2X,2HIS F12.5,2X, 
13HCPS) 
9808 FORMAT(1H1,3AX.43HWERE ARE THE EIGENVALUES AND EIGENVECTORS ///) 
RETURN 


END 
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FORTRAN DECK 


. § 
| CFULFL EXPANDS REDUCED FLEXIBILITY MAYRIX BY INSERTING 1 OR 2 ZERO 
c 


NIaQanrnr cs 


c 

1 

50 

96 

97 

96 

99 

100 


2 


10 


19 
18 


hy 


ROWS AND COLUMNS REPRESENTING ATTACH POINTS. 
THE ARGUMENTS ARE 
BCS) = REDUCED FLEXIBILITY MATRIX IN COMPACT FORM 
NxXC = ORDER OF RENUCER FLEX. MATRIX 
INPUT PATA RFQUIRED 
NR = NO. OF ATTACH POINTS (1 OR 2) 
NNE,NHO = KASS NUMBERS OF ATTACH POINTS 1 AND 2 RESP, 


SURROUTINE FULFL(BSNXC) 

DIMENSION B(1),N1(61275),0(50) 

DATA Q7/6HEXPAND/,OB/6HED FLE/.- O9/6HXIBILI/,OL0/6HTY MAT/.+911/6HRi 
1X / 

READ(5,1)NR,NNE,NWO 

FORMAT (918) 

MS=NXC4NR 

MMS=MSa(MS41)/2 

nO 50 121,MNS 

N(tT) 20.0 

JJJ=0 

KK=0 

Jus 

NO 100 Jsi,Ms 

IFC JeEQ.NNE.OR Je FO NWO)GO 10 99 

TesJei 

JIS LONXCe tJ 

KKK=J-1 

nO 98 JKFI,J3 

KKK=KKKel 

KK=KK+4 

TE CKKK EQ.NNE.LOR. KKK EO.NWO)GO 10 96 

60 10 97 

KK=KK+i 

NCKK ABC IK) 

CONT ENUE 

GO TO 100 

KK=KK*#MS-Jel 

JIJEI IJ e1 

CONTINUE 

WRITE(6,2) 

FORMAT(///86HU P PEAR 
JexiBpittidgty MA 
no 10 I=1,mS 
NS=(2eT+(]-1)#(2eHS-$))/2 
NE=(2eMS4+ (Tet aC 2eMS-1))/2 
WRITE(6,3)1,¢°00(5),J32NS,NE) 
FORMAT(/SHROW, 14/(9F14.5)) 
CONTINUE 

PUNCH 4,07,08,09,010,011 ¢ 
FORMAT (546) 

nO 20 1=1,MS. 

Tlsded 

IF CET.EQ-0) GO TO 18 

no 19 Jai,t) 

NU=(2e 14 ( Soi e(2eto- J) /20¢ U1 OCNS<1) 

C( J SDCNU) 

CONTINUE 

N§o(2efe(lo1)*#(2eMSel))/2 

NEx(2ehSe( fo] )e(2eMSef))/2 
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«cuneate aS LEIS SPTTEIS INSEE 


NGLE*EXPANDED Fee 


DO ON AV ®@ ON 


DO 16 JJ=NS,NE 
ecsen( as) 
16 J=U+1 
PUNCH 5(O(S),322,MS) 
> FURMAT(1P6€29,5) 
20 CONTINUE 
RETURN 
END 


$ tORTRAN BECK 

COTVID 

c N=NU. OF NORMAL RISPLACEMENTS ° 

C M=NQ. OF ROTATIONAL U.O.f. 

ct NIPE-CONTAINS STIFFNESS (OR ASS) MATRIX 

(' MIPE-K12 (412) STORED 

c TTPF~K13 (M11) STORED 

Cc A- DUMMY STORAGE VECTOR.-LARGER OF (N#(R+1)-2 OR M*(M41)7/2) 


SURROUTINE BIVID (NSM NTPE,MTPE,JTPE.A-B) 
NIMENSTON AC').BC1) 
REWIND ITPt 
REWIND NFP 
REWIfD MTPE 
NMAX=IN@ (Nt 572 
MMAX=NO (MO 8 72 
NM=N+M 
TONT=0 
aU JG F=1-N 
TESNMe-fo) 
REAVCNTPE) (805),J021-71) 
Tw=JI-M 
nO fl J=t,tv 
TUNTZIONTS: 
Pa ACTCNT)=B8¢0J) 
TVI=TH+) 
JOND = 
No of Jztp',tl 
AJCNTSJCNT?! 
a RC JONT)=BC0) 
WRITECMTIPE) (8¢4),J=7.M) 
i) CONTINUE 
WRETECT TPE) (4059, 521, NMAX) 
REWIND MTPE- 
REWIND PTPE 
Turk 
TUN] =i 
ho bt [si,m 
TL=M-fONT : 
READCNTPE) (8¢8),J35)-1910) ! 
TCNIZSCNT+: 
AO oN J). EF 
Toc} Nel 
on A(T =R( SI) 
»N CONTINUE 
RETURN 
ENN 
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FORTRAN DECK 


or TE ad 


ZROMAK 
D 1S A DUMMY VECTOR WITH STORAGE N OR MH (LARGER) 
A TS A BUMMY VFUTOR WITH STORAGE WNeC(N+tid/?2 UR Me Mel )s,2 CLARGER- 
B TS A DUMMY VECTUR WETH STORAGE N OR HM (LARGER) 
© TS A NUMHY VECTOR WITH STORAGE WN OR M (LARGER) 
N=N1}. OF NORMAL DISPLACEMENTS 
M=NU. OF ROTATIONAL 0.0.F. 
NIPE CONTAINS K1) MATREX 
MIPE CONTAINS Kitz MATRIX 
1TPE SCRATCH TAPE 
KTPE STORES KI) 2 #K-2868(-1) 
A INETTIALLY CONTAINS K°2 INVERSE 
#@e RENUCED STIFFRESS MATRIX IS STORED ON TIPE | 
SUPROUTINE ZROMAK(A,8,C,0,N,MNIPE.MTPE,LTPE, KIPE)D 
NIMENSTON 4¢1),8061),0(01),001) 
NUMHLE PRECISION SUM, DP1,vP2 
CALL SYMINVOC AoM) 
REWIND MTPF 
REWIND ETPF 
REWIND NTPF 
REWIND KTPE 
NMAX=N@(Ne+) 7? 
MMAK=M#(Me* )/7 
AQ 18 KK=1.N 
RE AUV(MTPE) (ROT), 721.) 
IUNT=0 
No 1.008 PK=',M 
JJ= IK 
WK 7K 
nu - 0 Jad. m 
TCNITSLONT+) 
on CCIE ACTONT) 
HaJJo1 
JA=M 
T0=TK 
AO oh Jet, dg 
TEC US-EQeuU) RA TO SU 
GCCJVSACED) 
JA=JA-1 
PU=ED+JA 
sa CONTINUE 
SUM=0.N 99 
no +a J=),M 
APIT=B(J) 
NPP=C( JS) 
a | SUM=SUM+DP: ape 
NC JK = SUM 
jhin CONTINUE 
WRETE CITPR) (OOS) S51, M)d 
WRITE (KTPE) (NOS) JEL M) 
1a CONTINUE 
REWIND ITPI- 


4 
4 
4 
H 
q 
} 
| 
4 


AARAIBAN 


- 
ry 
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pees ee 


REWINQ MTPE 
REWIND NTP 
REWIND KTPE 
READ (NTPE) (ACJ). J=1,MHAX) 
TUNT=9 
ng o8 KkK=) Pa | 
READ CITPE) (NCJ), 5=1,6) 
KI=KK 
ha 7fl KJ=1.N 
REAHCMTPE) (UCU), J=1.M) 
KP=kJ 
TECKPLLTKE) G0 TO 78 
SUM=0.9D0 
HO etl KR=1.4 
NPT=DCKR) 
AP? SE (KR) 
tn SUMESUM +DP, epP2 
ICNI=JONT+ 
SM=SUM 
ACTCNT)SACTENT)-SM 
70 CUNT ENUF 
REWIND MTPE 
of CONTINUE 
REWIND NTPE 
REWIND MTPE 
REWIND ITPE 
WRETFCITPE) (ACT). 1=1,NMAX) 
REWIND ITPF- 
RETURN 
FND 
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FORTRAN DECK 


CZRIMAM 


2a 


hf 


poe 


if 


N=NO. ‘OF NORMAL NISPLACEMENTS 

M=NU. OF ROTATIONAL D.U-F. 

NIPE CONTAINS Mili MATRIX 

MIPE CONTAINS M1? MATRIX 

{TPE SCRATTH TAPE 

KTPE CONTAINS Ki?eKs26e0e(-1) 

REDUCED BASS MATRIX IS STORED ON ITPE 
SURROUTINE ZROMAM(A.B.C,D No MNTPE.HTPE, TIPE. KITPE) 
DIMENSTON AC1).8013-C0¢01 3.0013 

DOUBLE PRECISION SUM1,SUMZ2-DPL. DP. DPS 
NMASS=N 

REWIND MTPE 4 
REWIND LTPF 

REWIND NTPH- 

REWIND KIPE 

NMAX=Ne@(NG6 72 

nO rth KK=L.N 

READ(KTPE) (ACT), 121-8) 

TONT=4 

NO 1000 JK75. 

JSOSSEK 

JK=1K 

NO 7 JaJJ-H 

[CNISICNT+! 

GC J)=aCTONT) 

JJeUJS<1 

JA=M 

1D=1K 

NO oft Jet, du 

Te¢(JJ-EQ-0) GO TO 34 

C(SFACTD) 

JA=JA~1 

TU=1D+Ja 


* CONTENUE 


SUM) =4. D0 

NO oh J=toM 

NP1=B(J) 

AP? =C(J) 

SUM) =SUMT FIP: DP? 

NC SK) =SUM1 

CONT ENUE 

WRITECTTPE) (00S), J=1,4) 
ONT TINUE 

REWIND TTPE 

REWING MTPE 

REWIND NTPE 

REWIND KTPE 

READCNTPE) (ACS), S21. NMAX) 
NO ol) KK=EL.N 

READCHTPE) (RC J). Jet 6M) 


112 


CE nn ee ecaaiiememmemeeciesamianemateentaraifaecnet ant cana taaat tant Te ca oe 


4 


REAUCETPE) (9¢9),522.4) 
HO ¢0 XJF1.N 
REAVOCKTPE) (005), 521.) 
Sums cu. bo 
SUM =i. Of 
NQ o0 KRF1. 
NP1=BCKR) 
NPP=DCKR?) 
NP %=C (KR) 
SUM) =SUMT +P 1s hPs 
wh SUB SSUMZA +H Poa PS 
Sai =SuMi 
SM? =SUM2 
TECK Gee KK) MM=( 2K JOCK Kat 8 CP aNMASS~KK)D/2 
TECRS.GEeKK) ACMM)=A CMM) =SM2E4+SMZ 
TRQCRJ.LECKK) MMS(- ORK EC KJ-1) OCP ONMASS KI) 72 
IF CKU.LEKK) ACHM)=A(MM)-SH1 
8 CONTINE 
REWIND KTPF 
HY CONTINUE 
REWIND NTPE 
REWING KTP 
REWIND ITP 
REWIND KTPEF 
URPIECETPE) CAC), P51 ,NMAX) 
REWIND FE TPH 
RE TURN 
ERD 


’ 
° 


b 
C 
C 
¢ 
c 


. 


C 
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FORTRAN DFCK : 


SYMINY 


tuft 


Vea 
aa 


Jon 


aun 
4¢e9 


500 


Aud 


Tun 
oun 
Taun 


A #8 THE UPPER FRIANGLE OF THE SYMMETRIC MATRIX fo RE 
FLFEMFEFNTS ARE STORED ROWWISE. 
N = ORNER OF HATRIX 
PROGRAM INVERTS IN PLACE. 
SURROUTINE SYMINV(A,N? 
NEMENSTON AC1) 
NMAX=NA(N0!) 77 
ACTI=SQRTCACT)) 

No itu JJ=v, 
ACEIDSACTJI/AG1) 
A(1)=1.0/A01) 

Im1=3 

IJ=N 

NO ifwsi f=-,N 

Piz t Jet 

tusif 

Nu 7fid Jet. 

JAT=eJ-1 

List 

LJeJ 

NO 170 LL. 1M) 

ALT JZACTJIAACL ID @ACLUD 
PStbonet 

LJ=LieJM! 

Tush Jet 
ACTI)SSORTCACET)) 

Jtsl 

JJz! 

NO bfu JF1-TML 
ACGIISACIUJI OAC UID) 

TE CU TMLD S96 4FU, dA 
JP1T=J+1 

JisJJ 

tLisJl 

NO 400 L=JZP1.11 
JL=JL+1 

LISLI+Neb +? 
ACJPIAAC JE OAC IJLIMACLI) 
ACJTS“ACJPDACED) 
Jt=J1leNe-J 

ARES Coe 

TE CT-N) 60. gun, OnL 
PP1=T+1 

TJ=l] 

nO shu J=IPI.N 

[JslJel 
ACTIDSACTJIZACTD) 

ACTED) =9.07ACTH) 

IM1=] 

T)=i 
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wee ee ee lll scemethemneratenn:tuemeenemenemensaentemaaneienetal 
ene SNR LTRS RT Ne 


INVERTED. 


AM LN NM. ANS. 


mM... MY 


-2* « 


M7 oR Jol 


oe MH oe wr WH 


170 


170 


1408 
eAhat 


00 chu THl,n 

SJ=1] 

tJ=1] 

DQ 140 J=t,Nn 

ACT JIHACE SI OAC IE) 
jP1=Je1 

TF CUPL=N)LIue,19GG07406 
tl=!J 

JLeJJ 

NO 1200 L=APION 

TL=ite1 

JL=sJb¢1 
ACTIJIZACTUIOAL EL) #ACIL) 
JJzJUb +) 

TJs} Jet 

ttetd 

RETURN 

END 
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SON Batt 


fORTRAN DECK iS 
TE EOMAT 
THIS SUBROUTINE FINDS THE EIGENVALUES AND EIGENVECTORS § OR 
SYMMETRIC MASS AND STIFFNESS MATRICES. 
THE ARGUMENTS ARE-- 
N- ORDER UF MATRICES. 
A~ DUMMY VECTOR WITH DIMENSION IN MAIN PROGRAM OF NeUNG1) 72 
VALU- STORAGE FOR FIGENVALUFS. MUST BE DIMENSIONED IN THE MAIN 
PROGRAM AS A VECTOR OF LENGTH WEIG. 
TEMP. H,C-U.- DUMMY VECTORS WITH DIMENSION OF N IN MAIM PROGRAM. 
F- DUMMY ARRAY WITH DEMENSIONS OF (N-3) IN MAIN PROGRAM. 
TOUM- DUMMY INTEGER VECTOR WITH DIMENSION OF N IN MAIN PROGRAK.) 
MIAPE- TAPE WHERE STIFFNESS MATRIX IS STURED IN COMPACT FORM. 
NIAPE- TAPE WHERE. MASS MATRIX IS STORES IN COMPACT FURM. 
JTAPE,TTAPF- SCRATCH TaPES. 
NE!G- WUMRER OF EITGRNVALUES DESTREL. 
NV¥EC- NUMRER OF EIGENVECTORS DESIRED. MUST HE EQUAL 10 OK LESS | 
THAN NEITG. 
THE WASS AND STIFFNESS MATRICES ARE STORFU IN COMPACT FURM AS 
VECTORS. ONLY THE UPPER TRIANGLE OF THESE MATHICES( UY RUWS) 15 
STORFD. 
SUBROUTINE EIGMATON, A, VALU. TEMP,8,C,.D,F+ TUM. MTAPE.NTAPE» JTAPE. 
LETAPE,NETG.NVEC) 
RIMENSTON A401). TEMPC1) VALUCED -BC2d-CC UT) DOLIZECNe 5), EDUMOD) 
DOUBLE PREC IESTON SUM. SUMi 
INTEGER OUT 
NUT=6 
REWIND ITAPE 
REWIND JTAPE 
REWIND NTAPE 
REWIND MTAPE 
Mav eh 
NMAX=Ne (NE) ) 42 
Ce 2s # #@ ¢ # © *# &# He HH oe He He eH He He er Hh le UU UU Ue Ue Ue 
im STEP 1 
C REAL IN M 4Y ROWS IN COMPACTED FORK 
C 
¢: 


SRI ASS Cees Sew ae eave narra Ae Sepa eee eoae ee 


REPLALE M KY (LX TRANSPOSE, WHERE M=Le(L TRANSPOSE 
CAICULATE FIRST ROW 
REA (NTAPH) (ACT) P= t+NMAX) 
REWIND NTAPE 

% CONTINUE 
ACT)=ZSARTCAC?H)) 
BO iff |2?.N 

18 ACEVEACT I AGT) 

C CALCULATE ALL THE OTHER ROWS 

TNION 
WO sf) FS¢e.N 
TND=IND+) 
SusM=0 Ne 
Kizbe-) 
HO OD JJELe KI 
MJ=(M~ SSD OC Set S74] 
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- Aa tet SR MNT A DLO LOE AAA 2 


ol SUM=SUM+A( MJ) AC HS) 
ACTND}=0SQRTCACIND)-SuM) 
TE CIND.EQ.NMAX) GO TO :UR 
SUMIZAC END) ; 
Ki=}4} 
nO v9 J=K1-N 
TND=ENDed 
SUM=0.04 
Tick 
na of JJ=L- TE 
K= (MSI HC Ie? 
KE=A+] 
KJrKeJ 
A SUB=SUMFACK EY OACK I) 
ACTND)=(ACINA)-SUM)/SUML 
29 CONTINUE 
tui CONTETNUE 
Tet CONTINUE 
c CHEEK FOR SINGULAR MASS MATRIX 
NO itv J=).N 
Ki=(M-fpeCt-')%2e] 
TECACKE) £0.29.) GO TO 1690 
14? NONTINUE 
c THIS COMPLETES STEP } 
Ce © # @ # # 8 # e # Fw ee © FE H He 8 HF HF HF FH HF FH &F + Hee F He HF HF H H F 
c STLP ? 
C WRITE (LITRANSPOSE ON FAPE BY COLUMNS 
C PUT (CLYTRANSPOSE INTO TFHPURARY STURAGE (TEMP~--A VECTOR) 
t ANU THEN WRITE TEMP ON TAPE 
KTAPE=NTAPE 
tiv ENN=A 
Nu 14 JELLN 
NU os TF1.J 
IND= ENS +t 
M)1s(Mep eC }-))72¢J 
TEMPCINDI=ACMII) 
Tol RLONTENUE 
WRITECK TAPE) (YEMP( JJ). JJst, END) 
TNP=0 
$a) CONTINUE 
THIS COMPLETES STEP + 
* » &*© # © #© © © # @&© # &©@ @ & © &@ © @ @ © & &#@ © & he He HH H He HH HH HO 


c 
Cc 
¢ STeP 3$ 
¢ 
c 


tt 


CCLITRANSPOSE) TNVERSF REPLACES (LY TRANSPOSF IN CORFE 
REPLACEMENT 1S PUNE BY LAST COLUMN FIRST--WORKING UP THE CULUMN 
nO atu TFl.N : 
PNM=(T#( M4 <1) 1/22N 
4c" ACTIND)=te/ACIND) 
NA aVy Jae. 
IS(N+2 DoJ 
ho .9n LF¢.J5 
INNS ON+ Jt] ~ eC Se ED/ 2 
SUM=09.DH 
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Kise yJ-1+2 
NO 650 KEK, Ud 
TDK=ITND+K 
MRE(MHK DOCK 1/24 IS a 
49 SUM=SUM*ACTIK) #ACMK 
THD: IND JS j 
EDE=ENN-1+! 
490 ACTND)=-SUH*AC LUT) 
499 CONTINUE 
t END OF STEP » 
ok a ee ee eZ 
Cc STEP 4 
c U=((L) TRANSPOSE) INVERSE 
¢ WREVE t ON TAPF RY ROWS 
WRETECITAPE) (ACT), 7=2t+NMAX) 
c FINISHED WITH STEP 4 
Ce & © # @ &# @ #© @ & © 8 & & & © e 8 ee HH eH HH HE HE HF HE eH HH HB 
STFP > 
WRITE U ON TAPE RY COLUMNS STARTING WITH THE LAST rc OLUMN FEIRSI 
PUT U (LAST CULUMN FIRST) INTU TEMP AND THEN WRITt OWN TAPE 
INN=0 
NY bS> KFe-N 
JEN-Ke) 
AQ “Su Frayed 
TNN=INDel 
MPP=(Me PAC Pet /P+S 
TEMP CINDI=A(CME?) 
5.8 RONTINUE 
WRITECJTAPE) (TEMP( JJ) JJ=te IND) 
INn=0 
5945 CONTINUE 


ans 


C END OF STEP » 
en ial 
C «STEP 6 

C FORM KU 

GC RFAN K INTO CORF 

© READ U INTO CORF A COLUMN AT A TIME IN REVERSE ORUER 

C  RFRLACE K BY KU COLUMN BY COLUMN STARTING WITH THE LASI COLUMN 

C AND WORKING UP THE COLUMN 


READCMTAPE) (ACL), 1T=1,NMAX) 
REWIND JTAPE 
NY o90 JJ= oN 
JENtT-JJ 
READCJTAPE) CTEMPCTT)-Plste J) 
NO 59 li=',J 
f=J+1-T1 
SumM=f1. Du 
NO O50 K=1.] 
MK1= (MK) (K-19/2 28 
46% SUMZSUM+ACMKT )*TEMPCK) 
INNS (Me PdeCF-LI77+S 
160} -EQ.J5) G0 Tu nBh 


wt 

% 
é 
se 
4 
4 
4 

E 
3 


Ki=(M-1T)@C F-19072 
12-141 
WO 66 K=J.Jd 
KIK=K1+K 
Bot SUMZSUMFACK LK DS TEMP CK) 
650 CONTINUE 
ACIND)=SUH 
490 CONTINUE 
END OF STEP 6 
7% & & © © # FF £ &F FF FF &F HH FF RF HF HH He FH FH FH He HF 8 FF FH Ht 
STEP 7 
FORMC (LOI ENVERSF ) ok 
KU TS EN CORE 
RLAD IN L EOLUMN SY COLUMN AND CALCULATE CCL ENVERSE) #KU 
ROW BY ROW 
CALCULATE THE FIRST ROW 
REWIND NTAPE © 
REANCNTAPE) TEMPC 1) 
NO /ia [=t-N 
7i0 ACTISACEO/TEMPOL) 
c NOW CALCULAIF THE REST OF THE ROWS 
IND=N 
N09 /99 [=2.N 
REAL (NTAPE) (TEMP( JS). JJs151) 
NO 499 JFT-N 
TND=TNDel 
JJsl-t 
SUM=0.) 
NG sou KEL dd 
MKZ=OMH—K DC K-19 /765 
7on SUM=SUM+ TEMP CK) *A CHK? ) 
799 ACEND) = (CACIND)-SUM)/TEMPC 1) 
STEP 7 IS COMPLETE 
*# &¢ # #@ & # &© #* &€ # &@ 8 © © & & # &© & 8 F H # HH 8H He FH HR He HH HF HF H eH H OH 


¢ 
C 
C STEP 8 
Cc 
C 
Cc 


nan 


AINA BT 


DETERMINE EIGENVALUES AND EIGENVECTORS OF THE NEW MAFREX 
CHANGE THE STGN OF A IN ORDER TO URTAIN THE SMALLEST 
FIGENVALUE FIRST 
NO eh JFLeNMAX 

Aug ACT )=-ACT) 
CALL BIGMATCA, VALU. TEMP,8,0,D,E>-IDUMLN-NEIGSNVEC-MTAPE) 
c CHANGE VALU RACK 
NO o5u0 TEL. NEIG 
850 VALUCE) S-VALUCT) 
( STEP 8 IS COMPLETF 
C# 2 # © & # © © © # © & & &© & * He He He Hw % FH FH He eH eH HY & He HH F H DB * 
c STFP 4 
C CHANGE EIGENVECTORS BACK 
c REAN U INTO CORF RY ROWS 
C READ UNCHANGED FIGENVECTORS INTO CURE ONE AT A TIME 
c CHANGE AND PRINT EIGENVECTORS 
TH CNVEC. EQ.) GO TO PuttU 
WRETECOUT, burt) 
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REWIND ITAPE T 
READCITAPE) (ACI), 1=1,NMAX) 
REWIND ATAPE ; 
NE ¥Oy JJ=',NVEC { 
READCMTAPE) (TEMPCT).121,N) 
IND=O 
NQ Yiu FSien 
SUM=H. AD 
nN YOO JEN 
IND=IND+1 : , 
Ou SUM=SUM*ACIND) © TEMPCJ) j 
9if TEMP(})=SUM 
f NOQMALIZE THE EIGENVECTOR 
SUM=TEHP(1) 
NU 939 ETON { 
TF CABS(SUM)~ABS(TEMPCIEID)) 950.939.9349 
Qo SUMSTEMPTIET) 
939 CONTENUE [ 
TF ESUM) 94°-,947,94N 
940 CONTINUE 
Nu v4. Tl= oN . 
TEMPC EE) STEMPCYTIT)/ SUM 
94) CONTENUE 
947 CONTINUE 
999 WRITECOUT, 2b tu) SUs VALI SJ), (TEMPCE), 1=1,N) 
c STEP 9 IS COMPLETE 
Ce #*# eee Oe ee eee eee a 
GO 10 PuoY i, a l 
4itut FURMAT (1H ",'OH EIGENVECTOR MUMBER 15/3°X.,17H CORRESPONDING JD 
TIPELSSI/CLH TPAEI 227 3) : 
41 FORMATCIHL. SSX. 43HHERE ARE THE CIGENVALUFS ANU LIGFNVECTORS fil, 
400? FORPAT(IVHI].SYX.27HTHE MASS MATRIX IS SENGULAR 4/7) | 
19ye" WRETFCOUT, -& 7) 
200" RETURN 
ENR 


 ecinateiendshaamaitinetin mae somes sot tuk. ones coaameien a 


wot. Kebd wtih ote. 


$ FORTRAN DECK 

CHIGHAT 

C PROG.AUTHORS M.ELSON AND R.~E.FUNDERLIC-CENTRAL DATA PROCESSING» 4.14.05 b 
SUBROUTINE BIGHAT(As VALU, VALL,UPPERU, DIAG. V,1T, INTER, NN NEIGSNVFL, 


* 
Po, Ee) ee PSE sk RE, NPN Dee ae UR GS ee 


IMTAPE) 
NIMENSTON 4¢7),VALU(1), VALE C1). UPPFRDCI) BIAGLT).VCL)>TONNDO)- 
LENTERCI) 
REWIND KTAPE 
NZ=t H 
N=NN n 4 
IF (NSLF.2)60 TO 49 , 
NPI=Ne1 kK 3 
NMT=N-1 RO 3 
NM?=N-2 ts : 
NI? PL SNe? 4: R ; 
TX=t t 
NO 30 T=i,NM? rf) 
SIGMA? =. H 
[PV] +4 ] 
NO 5 JeIPL.n 2 
PJU=}XeJ K 
1 SIRMAC=SIGMA FACT J) 887 kK 
SIGMA=SORT(SIGKA?) H 
PisiXe] ts 
DIALCT)=SACTID i 
TIPtsyX+de! H 
UPPERD(E)=-SIGN(SIGHA,ACTIP1)) 7 
T(1,7)=S9TGMA? 
Yt CABSCSIGMA).GI.ABSCACTIP1)2)66 FO 2? 
HPPERDCJT=ACTIPL) H 
ACTIPL) =f. H 
60 10 18 ) 
> ACTEPL=SORT (C1. +anSCACTIPL))/SIGMA) b 
SQTOAM=-SIEBN(SIGMASACTIP1),UPPERDCI)) b 
PP?=]+? kB 
nO » J=TPe.N Kk 
Tu=IXed 8 
S ACTJZACES)/SQTGAM ) 
ARP eTa(2@nN- =) 977 iH 
IX= UKE is 
TIX=JK1 H 
N0 % JSTPL.N 3} 
VALL(J)=0. 
JK=IKL 6d K 
00 4 KSIP1.J ts 
[K=1X+K ' H 
VALE CSIEVALL CJD OAC IK DOA CIK) 
A SK= JK eNe-K 8 
TF (JeEO.NIGO TO 6 Ki 
CALL LOOPL(C U2. NP1>VALL (JS) -ACIX)-ACEX)) 
h JX UK+Ned Kk 
& NELGAMZI. b 
NO . JEIPL.N R 


rh 


Wy 


i% 


14 
19 
tfé 
L7 


wR 
jo 


en 
a1 
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